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EXPERIMENTAL STUDIES OF THE VAPORIZATION 
OF DROPLETS IN HEATED AIR AT HIGH PRESSURES 
Clyde William Savery 

Under the Supervision of Associate Professor Gary L. Borman 

ABSTRACT 

Measured histories of vaporizing n-heptane and Freon- 
13 droplets suspended in a heated air stream are reported. 
The range of air conditions was 1.5 < P < 100 atm and 
100 < T ra < 300° F. Comparisons are made with film theory 
calculations corrected for the effects of total pressure 
on thermodynamic properties. Agreement was satisfactory 
for liquid steady-state temperatures and mass transfer 
rates except for Freon-13 where significant differences 
were found at air pressures above 1.25 P_. Freon-13 droplet 
histories were found to be completely unsteady for pres- 
sures above 60 atm. The absence of steady-state conditions 
was found to correlate with the critical mixing point for 
the Freon-13-air system. 



SUMMARY 


The problem of droplet vaporization at ambient gas 
pressures in the near and supercritical pressure range of 
the vaporizing liquid was investigated by a technique where- 
in drops were suspended on a temperature probe and observed 
by means of shadowgraph movies. The forced flow of heated 
air was upward and uniform during the droplet lifetime. 

N-heptane drop vaporization data corresponding to the 
range of conditions, air temperature 100 - 300 ° P, air pres- 
sure 1.5-100 atm abs, drop Reynolds number 150-600, and 
initial drop diameter 1275-2025 microns were obtained. 

These data constituted the low ambient gas temperature 
(T oo /T c < 0.8) and high pressure (0 < P/P Q < 3*7) range for 
n-heptane which has a critical temperature of 312° P and 
a critical pressure of 27 atm abs. 

Chlorotrifluoromethane (Freon 13) drop vaporization 
data corresponding to the range of conditions, air tempera- 
ture 100-300° P, air pressure 29-68 atm abs, drop Reynolds 
number 250-450, and initial drop diameter 980-1200 microns, 
were acquired. The Fr-13, which has a critical temperature 
of 84° P and a critical pressure of 38.2 atm abs, results 
are in the range of moderate gas temperature (1 < T^/T < 1. 

and high pressure (0.75 < P/P < 1.75). 

c 



Vaporization rates and drop plateau temperatures pre- 
dicted by the application of existing empirical heat and 
mass transfer correlations were compared with measured 
values for individual droplet histories. The calculations 
accounted for high pressure vapor-liquid phase equilibrium, 
enthalpy of vaporization at high pressures, drop surface 
regression, and high mass transfer rate disruption of pre- 
dicted heat and mass transfer rates. Calculations with 
these high pressure corrections produced agreement with 
the measured values to within the estimated comparison 
uncertainty range of + 13° P for plateau temperature and 
+ 25$ on mass transfer rate for all data except that in 
the moderate ambient gas temperature range and reduced 
pressures above 1.25 where the calculations underestimated 
mass transfer rates by more than 25$ ■ Calculations with 
low pressure assumptions gave good predictions of drop 
plateau temperatures over the entire range but underesti- 
mated mass transfer rates considerably at gas pressures 
above 5 atm abs. 

A departure from the commonly observed drop tempera- 
ture history phenomenon wherein the temperature rapidly 
approaches a plateau value and remains there until vapori- 
zation is practically complete was observed at moderate 
gas temperatures and reduced pressures greater than 1.5. 
This new vaporization regime which is identified by a 



steadily increasing drop temperature was successfully 
correlated with the drop interface gas mixture critical 
state in the range of reduced gas pressure of 1.5-1. 8 
with respect to the critical pressure of Freon-13. 
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INTRODUCTION 

The importance, of droplet vaporization to the combus- 
tion process in diesel, jet and rocket engines has been 
realized for many years. Early studies of this phenomenon, 
both in chemical engineering and combustion, were confined 
to low pressure conditions. Several recent articles review 
the data accumulated since World War II for heat and mass 
transfer from single drops (19, 25, 50).* 

In the last decade improved performance in combustion 
devices has led to ever higher operating pressure levels 
which are now approaching the critical pressure range of 
fuels. In addition, oscillatory combustion in liquid- 
fueled rocket engines creates large amplitude variations 
in the instantaneous chamber pressure causing droplets to 
experience much higher - than normal design pressures. 

The ability to predict both steady state and transient 
performance of engines is dependent in part on the success- 
ful application of a droplet vaporization model in the high 
pressure region. The application of empirical heat and mass 
transfer coefficients based upon boundary layer average 
properties led to a film theory vaporization model which has 
been successfully .correlated with data at atmospheric pressure 
(6, 44). Thus, strong motivation exists for adapting the 
film theory model for the high pressure region and testing 
it with the results, of high pressure vaporization experiments. 

^Numbers in parenthesis refer to references listed alpha- 
betically by first author. 
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At atmospheric pressure both experiment and theory 
show that a drop which is suddenly exposed to a high tempera- 
ture inert gas environment will first show a gradual rise 
in droplet temperature followed by a period of constant 
droplet temperature. The rising temperature portion of 
the history is referred to as the unsteady or "heating up" 
period and the constant temperature portion as the steady 
state period. The steady state temperature will be referred 
to here as the plateau temperature or "wet bulb" tempera- 
ture. During the unsteady state the droplet temperature 
rise causes the mole fraction of droplet vapor at the 
droplet surface to increase, thus providing a larger driving 
potential for mass transfer. At the same time the decreasing 
temperature difference between droplet and ambient gas 
causes a decrease in the driving potential for heat transfer. 
During the unsteady period some of the energy reaching 
the droplet surface goes into vaporizing the liquid, and the 
excess goes into sensible heating of the droplet liquid. 

From previous studies it has been shown that because of 
internal circulation within the drop, the droplet tempera- 
ture is essentially uniform being then only a function of 
time. At steady state all of the heat transfer reaching 
the surface goes into mass vaporization and the droplet 
temperature remains constant. The ratio of heat to mass 
transfer coefficients is essentially independent of rela- 
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tive gas velocity and droplet size. Thus, for a given 
liquid and given ambient gas, the plateau temperature 
depends only on ambient gas temperature and. pressure. 

Turning to the effects of high ambient pressure on the 
vaporization process, one must first recognize the existence 
of critical points for binary mixtures. At high pressures 
the ambient gas dissolves in the drop forming a mixture in 
the liquid phase. A critical point for a two component mix- 
ture of a given composition is the single point on the line 
bounding the two-phase region of a pressure-temperature 
diagram where the phases become indistinguishable (23). 

The locus of these critical points for all compositions is 
called the critical mixing line. If one assumes thermo- 
dynamic equilibrium at the droplet surface, it is intuitive 
that at sufficiently high ambient temperatures and pressures 
the mixture at the surface may reach or exceed the critical 
mixing point state. Such states usually occur for ambient 
pressures above that of the pure liquid critical pressure 
and mixture temperatures well below the critical temperature 
of the pure liquid. 

For droplet histories where the mixture critical point 
is approached, the thermodynamic and transport properties 
become increasingly dependent on pressure. In particular, 
the mole fraction of droplet vapor at the droplet surface 
and the enthalpy of vaporization must be calculated to inclucje 
the influence of total pressure. The enthalpy of vaporization 



4 


(latent heat) vanishes at the mixture critical point. Thus, 
if the droplet interface conditions reach the mixture 
critical point during the unsteady portion of the history, 
all subsequent heat transfer to the drop will result in 
sensible droplet heating. For such conditions no plateau 
temperature will exist, and the droplet temperature history 
will be completely unsteady. 

Recently two investigations of high pressure drop- 
let vaporization have been undertaken at the University of 
Wisconsin under NASA sponsorship. In the first, Manrique 
performed detailed calculations of the vaporization of 
carbon dioxide droplets in a high pressure nitrogen atmos- 
phere (35). Both steady state and transient vaporization 
models were investigated as well as the case of pressure 
fluctuation. In the second investigation, which is reported 
here, a high pressure chamber which can achieve air tempera- 
tures up to 300° F and pressures up to 100 atm was built to 
study suspended droplets. An upward flowing stream of heated, 
high pressure air vaporizes the droplet. The experiment 
gives droplet temperature histories and two-dimensional 
size histories via motion picture records. 

Water, n-heptane, and Freon-13 (Fr-13) data were 
obtained. The apparatus was tested with water drops. The 
n-heptane experimentation provided data in the region of 
supercritical pressure but at temperatures below the near 
critical range and below the ignition point. Fr-13, 
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chlorotrifluoromethane , was selected because Its low criti- 
cal temperature, 8M° P, gave promise of achieving a critical 
state during vaporization in the rig, and its high thermal 
stability with air avoided consideration of chemical reac- 
tions or combustion. In addition, Fr-13 is relatively 
non-toxic and has a moderate critical pressure of 38.2 atm. 

The high pressure data are compared with calculations 
based on a film theory that has been modified for applica- 
tion in the high pressure region. The transition from 
quasi-steady to unsteady drop temperature histories is 
examined, and its relation to the critical state of the 
gas mixture in the boundary layer is demonstrated. 

The combustion vaporization literature is reviewed 
in Chapter I. In Chapter II the film theory which was 
developed with low pressure assumptions (17) is presented, 
discussed and modified for high pressure application. 

Chapter III contains a description of the experimental 
apparatus and techniques. 

The experimental problems encountered and the data 
reduction techniques employed are discussed in Chapter IV. 
The experimental results are presented in Chapter V where 
they are compared with calculated values. Finally, in 
Chapter VI the conclusions of the study and recommendations 
for future work in this research area are given. 
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I. REVIEW OF THE LITERATURE 

Droplet vaporization has been an active area of com- 
bustion research since the early 1950’ s. The earlier work 
which is briefly summarized below was concentrated in the 
low pressure region. 

As time has progressed, attention has been increasingly 
focused on the high pressure region. In 1959 Spalding (53) 
proposed a droplet burning rate model for high pressures. 
Recent experimental results of burning droplets first gave 
evidence of a new mechanism of vaporization at gas pressures 
exceeding the critical pressure of the vaporizing liquid 
< 1 8 ) . Although the mechanism was recognized by the experi- 
menters and correctly associated with the critical region 
of the vaporizing media, a firm physical basis has only 
very recently been tendered. These interesting develop- 
ments are reviewed in this chapter. 

Experimental and Theoretical Results at Atmospheric Pressure 

A wealth of data on the vaporization of suspended drop- 
lets at atmospheric pressure was obtained by Priem and his 
co-workers (44). The droplet temperatures and size histories 
data covered a range of conditions : air temperatures of 

680-1430° R, air velocity of 27-228 in. /sec, initial drop 
diameter of 490-2180 microns, and 10 paraffin fuels. They 
compared their data with histories calculated with quasi- 
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steady low pressure, film theory. 

The film theory used by Prlem was first developed by 
El-Wakil, Uyehara and Myers (17). The Ranz-Marshall corre- 
lations for heat and mass transfer provided the empirical 
heat and mass transfer coefficients. The coefficients as 
applied include corrections for unidirectional molal diffu- 
sion and the effect of mass transfer on heat transfer. Ideal 
properties such as the pure liquid latent heat of vaporiza- 
tion and vapor pressure relationship were assumed. 

In applying the low pressure film theory model to his 
data, Priem found that temperature histories agreed to with- 
in 10° P and mass transfer rates within 20$. The agreement 
was considered to be within experimental accuracy. 

Experimental Results at High Pressures 

The earliest vaporization data in forced flow at pres- 
sures above atmospheric of which the author knows is that 
of Ingebo (26). Ingebo’s data included four pure liquids 
vaporizing in air at pressures of 0.2-2 atm abs and methanol 
vaporizing in air, argon, helium, and carbon dioxide at 
atmospheric pressure. A steady state porous sphere arrange- 
ment was utilized in the study. He found that the film 
temperature difference and mass vaporization rate both de- 
creased with increasing pressure, other conditions remaining 
constant . 

Hall and Diederichsen (21) using an apparatus , ’-which 
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was built for the study of suspended burning drops, obtained 
evaporation lifetime measurements of water drops in air 
heated by a gas flame or a carbon arc igniter. The pres- 
sure range of 1-20 atm was covered. They found that the 
rate of vaporization depended on the second power of drop 
radius and that the evaporation lifetime of a given diameter 
drop decreased with increased pressure. 

Petrazhitskiy (42) reported the results of an exten- 
sive experimental investigation of droplet vaporization at 
high pressures. He obtained liquid plateau temperatures 

j / t >-2 \ 

and vaporization constant measurements, — for water, 

9 6 % ethyl alcohol and kerosene drops suspended and vapori- 
zing in flowing, heated air. The air pressure was varied 
from 11-41 atm in each test except for water where the 
range was 11-60 atm. A Reynolds number range of 100-500 
was covered in the tests. His quasi-steady results with 
measured flow conditions were given with neither compari- 
sons to theory nor developed correlations. 

In their combustion tests of aniline droplets at high 
pressures, Brzustowski and Natarajan (7) performed some 
vaporization experiments. They vaporized supported aniline 


droplets in a mixture, : of 


nitrogen and 10# oxygen. The 


drops were located injat-closed chamber and in one series 
of tests were heated by carbon electrodes and in the other 
by the combustion of a solidr,propellant>- strand. The maxi- 
mum temperature near the drop and the evaporation constant 



9 


were determined. The test chamber pressure was varied 
from 100 to 800 psig. They concluded that the evaporation 
mechanism of aniline into a gas above its critical pres- 
sure but below its critical temperature is similar to its 
mechanism at low pressure. 

A fifth series of high pressure vaporization experi- 
ments was performed by Torda and Matlosz (55). They ob- 
tained temperature and size histories of n-pentane droplets 
injected into and supported in a heated and pressurized 
test chamber containing nitrogen. Data covering a range 
of pressure of 200-1400 psi at a fixed gas environment 
temperature of 850° R were reported. Presumably, the mode 
of heat and mass transfer was a combination of free convec- 
tion and simple diffusion. Their attempts at simulating 
the experimental results with calculations were not particu- 
larly successful. 

Borman et al. (6) reported the results of suspended 
hexane, decane, and hexadecane drops vaporizing in air at 
various temperatures and pressures from 1 to 5 atm abs with 
measured flow conditions. The experimental drop tempera- 
tures and mass transfer rates compared favorably with com- 
puted results over the rang e of temperatures and pressures 
investigated. 

Of all the high pressure vaporization experiments 
summarized above, only the work of Brzustowski and Natarajan 
and that of Torda and Matlosz extended into the supercritical 
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pressure range. Of these two experiments., apparently only 
Torda and Matlosz .observed steadily increasing droplet tem- 
peratures which did not attain the plateau which is charac- 
teristic of low pressure drop temperature histories. All 
six of the vaporization experiments summarized above 
utilized supported drops in a 1 £ gravity field. 

NASA Lewis supported work is in progress at Illinois 
Institute of Technology and University of Wisconsin on free- 
falling drops vaporizing in high pressure gases. The prin- 
cipal investigators of these projects are Professor T. P. 

Torda at IIT and Professors P. S. Myers, 0. A. Uyehara 
and G. L. Borman at Wisconsin. 

A few studies of burning drops at high pressures have 
been reported. In addition to their preliminary vaporiza- 
tion experiments previously discussed, Brzustowski and 
Natarajan (7) carried out burning experiments with suspended 
aniline drops at pressures up to slightly above the critical 
pressure of aniline. They inferred from their motion picture 
data that a different regime of burning occurred in the super- 
critical pressure range. Hall and Diederichsen (21) ob- 
tained burning lifetimes of suspended kerosene, tetralin, 
decane, furfuryl alcohol and amyl acetate drops at elevated 
pressures but not into the near or supercritical pressure 
region. 

Probably the most interesting high pressure droplet 
burning experiments, at least from the standpoint of vapori- 
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zation mechanisms, are those of Paeth, Dominicis and Olson 
(18). They measured the temperature and size history of 
suspended burning droplets In a zero £ pressurized test 
chamber. In a series of tests ofn-decane burning at in- 
creasing air pressures, the temperature histories showed 
a decreasing tendency toward a plateau attainment, and at 
a reduced pressure of about 1.6, only a slight inflection 
point in the temperature trace occurred. At the highest 
reduced pressure reported, 3*5, the droplet temperature 
increased steadily. They also noted the tendency for com- 
bustion lifetime to increase with pressure at pressures 
somewhat above the critical pressure. 

Theoretical Results at High Pressures 

Spalding (53) raised doubts about the application of 
the low pressure quasi-steady burning theory to predict 
droplet burning times at high pressures. In the high pres- 
sure region he abandoned the quasi-steady theory wherein 
the burning rate was vaporization rate limited. Spalding's 
theory considered the fuel to be a point source and utili- 
zed a first order chemical reaction as the outer boundary 
condition. He obtained a penetration theory result for 
the transient constant property case. As an initial 
condition, he assumed that all of the fuel was introduced 
instantaneously. His theory which is inapplicable to 
purely vaporization phenomena because of the outer boundary 
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condition Is Interesting nevertheless since it avoids the 
problems of transition to the critical state and the corres- 
ponding vanishing of the enthalpy of vaporization. The 
point source assumption has recently been relaxed by Rosner 
( 49 ). 

Wieber (57) investigated the attainment of droplet 
critical state conditions in rocket chambers. The results 
of his quasi-steady film theory calculations indicated 
that oxygen and heptane droplets would approach their cri- 
tical temperature at reduced pressures of approximately 
one and two, respectively. Ambient conditions of 5000° R 
gas temperature, 800 ft /sec final gas velocity, initial 
conditions of 10-1000 micron drop size and 80 ft/sec injec- 
tion velocity were considered. Furthermore, the trends 
of his results showed that the fractional droplet mass 
vaporized before attainment of critical temperature decreased 
with increasing pressure. 

More recently, Manrique and Borman (36) investigated 
high pressure vaporization. They performed detailed numeri- 
cal calculations for the steady-state vaporization of carbon 
dioxide drops in a nitrogen environment. For the case of 
a spherically symmetric boundary layer (i.e., no free or 
forced convection heat and mass transfer), they included 
the effects of thermodynamic non-idealities, property varia- 
tion through the boundary layer and high pressure liquid- 
vapor equilibrium. 
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The most significant result of Manriq.ue ' s and Borman's 
analysis is that at sufficiently high pressures they found 
steady state conditions cannot be obtained. Prior to their 
work it had been realized that if the liquid component 
could reach its critical temperature, a different regime 
of vaporization would be expected. However, they focused 
attention on the state of the gas mixture at the liquid- 
vapor interface. They pointed out the fact that this mix- 
ture could reach its critical state at a temperature below 
the pure liquid critical temperature but above the pure 
liquid critical pressure. This criterion determines whether 
a new regime is obtained. It would be expected that pro- 
bably a third regime of droplet vaporization would be reached 
when the droplet temperature reached Its critical tempera- 
ture and a puff of supercritical fluid resulted. 

Manrique (35) has extended his work to the unsteady 
case and has made comparisons with quasi-steady film theory 
models. Generally, he has found that the simpler theories 
give a fair prediction of vaporization lifetimes but give 
fairly large errors in predicting temperature and mass his- 
tories. Under some conditions of pressure and temperature, 
he found that film theory gives physically unrealistic 
steady state results. In the high ambient density region 
(high pressure and low temperature), film theory corrected 
for high pressure produced better results, than the low pres- 
sure model. However, the low pressure model gave better 
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predictions than the high pressure film theory model in the 
high temperature, high pressure region although neither 
model gave close agreement with the exact calculations in 
this region. 

In view of the recent work, two regions of droplet 
vaporization behavior were expected: 1) subcritical pres- 

sure range exhibiting a quasi-steady drop temperature 
history; 2) supercritical pressure range exhibiting a 
completely unsteady drop temperature history and charac- 
terized by a supercritical gas mixture in the boundary 
layer. The droplet temperature history could be expected 
to exhibit the classical plateau only in the low pressure 
regime. Departures from this form would indicate a transi- 
tion to the second regime. 

Manrique's calculations suggest that film theory modi- 
fied for high pressures may be expected to give approximate 
results at pressures in the near critical region. In addi- 
tion, a correlation between transition from the first to 
the second regime and the gas mixture critical point should 
be obtainable. These predictions and modifications to the 
low pressure film theory are discussed in the following 
chapter. 
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II. DISCUSSION OP THE THEORY 

In this chapter a physical picture of droplet vapori- 
zation Is presented. This is followed by the development 
of the film theory equations. The topics of gas mixture 
critical state calculations, property non-idealities due 
to high pressure effects and spatial variations through the 
boundary layer are discussed in relation to the film theory 
corrections required at near critical conditions. 

Physical Description of Droplet Vaporization 

The process of droplet vaporization can be idealized 
as a boundary layer problem with simultaneous heat, mass 
and momentum transfer between a stationary sphere and a 
steady flow of gas. 

Boundary layer development times equal to droplet 
radius divided by gas velocity were found to be short 
compared to droplet lifetimes, thus allowing the steady 
state idealization (9). 

Contrary to the usual boundary condition in rigid 
sphere cases, i.e., zero velocity in the zimuthal direc- 
tion at the sphere-boundary layer interface, a finite velo- 
city results from droplet internal circulation effects. 
Winnikow and Chao (58) in their studies of liquid-liquid 
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systems have correlated their observed reduction in measured 
drag coefficients with an internal circulation paramfeter, b, 
which involves the viscosit ieseof the dispersed (d) and 
continuous (c) phases. 



They also observed an increase in the angle of the boundary 
layer separation point and its dependence on their circula- 
tion parameter. Predictions of the effects of internal cir- 
culation upon drag coefficients in any systems other than 
highly purified fluids are, however, very questionable 
because of the large dependence of interfacial behavior on 
small impurities. 

When the drop Weber number, the ratio of drag to sur- 
face tension force, reaches a critical value, droplet oscil- 
lations begin. If no stabilizing trends ,such as radius re- 
duction due to vaporization or relative velocity reduction 
due to droplet momentum gain are established, droplet break- 
up in the case of unsupported drops or loss of droplet sup- 
port in the case of supported drops occurs after a time lag. 
Results of critical Weber number experiments, vary, but a 

typical value is 4 where the Weber number is based upon 

/ 

diameter (58). This criterion, although providing insight, 
does not strictly apply to supported droplet behavior. In 



addition to the drag and circumferential surface tension 
force on a free drop, a supported drop has a surface ten- 
sion force at the support point. The latter force accomo- 
dates itself to balance the weight' and drag force and can, 
in fact, have a resultant force either in the direction 
of the flow or in the opposite sense. Increases in droplet 
temperature and pressure will increase droplet Weber number 
through the effect of decreasing surface tension thereby 
increasing the possibility of oscillation. 

The conditions of the experiment and, in general, those 
existing in sprays, place the flow in the laminar regime. 

It would be possible to calculate the heat and mass transfer 
rates by solving the boundary layer equations except for the 
effects of flow separation. For average values of heat and 
mass transfer coefficients, one can resort to semi-empirical 
correlations such as those of Ranz and Marshall (46). The 
data synthesized with the Ranz-Marshall correlation were 
obtained under the conditions of low mass transfer rates, 
and the heat transfer results were independent of mass 
transfer. These correlations can be utilized with appropri- 
ate correction (4) for simultaneous heat and mass transfer 
at high mass transfer rates. 

Normally, only, two parameters, Reynolds number and 
Schmidt number, are used to correlate the dimensionless 
mass transfer coefficient, Nu A g(Re,Sc) . However, at high 
mass transfer rates in fluid-fluid systems, one would expect 
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a more complicated Sherwood number dependence: Nu nr , = 

f(Re, Sc, B, b). The B parameter accounts for the effects 
of high mass transfer rates (1). The effects of internal 
circulation would be accounted for by the b parameter 
given by Chao and repeated in Eq. (2.1). 

Film Theory. 

Film theory applies dimensionless transfer coefficients 
based on average film properties to compute heat and mass, 
transfer rates. In this section the film theory differen- 
tial equations for droplet temperature, T^, and mass, m, 
as a function of time, 0, are derived .under the assumptions 
of a quasi-steady model. Then the steady state result is 
developed as a special case. Physically, the steady state 
result corresponds to a simulated sphere experiment where 
the vaporization rate is balanced by the liquid feed rate. 

As a preliminary the definition of the mass transfer 
coefficient, k , is given by Eq. (2.2). Its definition is 

A 

based upon the molar flux of species A with respect to the 
molar average velocity. This convention, which is in accor- 
dance with Bird, Stewart and Lightfoot (4), has the advan- 
tage that k has a simple dependence upon concentration level 
and mass transfer rate.. In other words, k x gives the mass 
transfer coefficient for mass transfer alone; the bulk- 
flow term is subtracted. 

. ■ N AO - W N Ao + *Bo> 

x y Ao " y A» 


( 2 . 2 ) 
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N Ao , Ng o and y Ao are the molar flux of A, the molar flux 
of B and the mole fraction of A at the droplet surface and 
y A<jo Is the mole fraction of A In the free stream. If the 
ratio of fluxes, N R , is defined as 


N R " N Bo /N Ao> 


Eq, (2.2) becomes 


(2.3) 



W 1 - y Ao (1 * N R )] 

y Ao - y A~ 


(2.4) 


The superscript black dot (•) denotes that the transfer 
coefficients depend on mass transfer rate. The film theory 
summarized in paragraph 21.5 of Ref. (4)* was selected to 
relate the coefficients h* and k to the transfer coeffi- 
cients, h and k , which apply at low mass transfer rates. 

A 

The film theory corrections were obtained by assuming 
laminar flow in plane geometry. The fluxes, in terms of 
the boundary layer thicknesses, resulted from integrating 
the equations of change assuming constant properties. Then 
the boundary layer thicknesses were evaluated in terms of 
the uncorrected transfer coefficients. There are two major 
assumptions inherent In the derivation (4) which may have 
Implications for application to high pressure droplet eva- 
poration. The first and probably most serious is the 

#This <f film theory" uses a simplified unidirectional trans- 
port model to predict the variation of the transfer coeffi- 
cients with mass transfer rate, and it should not be confused 
with the use of the term in the present work. 
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assumption of constant properties 1 across the boundary layer, 
particularly constant heat capacity and the neglect of heat 
of mixing between the vapor and continuous gas phase. The 
second questionable assumption is that the film thickness 
is unaffected by high mass transfer rates. 

The results of these film theory correction are given 
by 

h' - h — - — (2.5) 

e - 1 



mu + r ab ) 


R 


AB 


where 


z _ N Ao (C pA * *R C pB> 


( 2 . 6 ) 


(2.7) 


- y Ao ~ y A<” 

AB " 1 

1 + N R “ y Ao 


( 2 . 8 ) 


and C pA is the vapor specific heat of species A. The para- 
meter z is equivalent to the "rate factor" for heat transfer 
and R AB to the mass transfer "flux ratio" cited in Ref. (4). 

With these preliminaries in mind, the quasi-steady 
energy equations can be written for a vaporizing drop. There 
are a number of assumptions implicit in this development, 

1. The drop temperature at any moment is uniform. 
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2. Viscous dissipation is negligible. 

3. Work done by pressure or gravitational forces 
is negligible. 

4. The energy transferred to the drop by radiation 
is negligible. 

5. Thermo- diffusion and diffusion-thermo effects are 
negligible. 

6. Steady state heat and mass transfer correlations 
are applicable. 


The following word equation corresponds to the terms in 
the energy balance, Eq. (2.9) 


- 


- m 


■> “ 

Energy 

stored 

s 

Energy transferred 
to drop by 


Energy carried 
away from drop 

in drop 


convection 


by mass transfer 


7T 


dT 

p L C pL d6 


L _ 


= it D h 


(Too - V - 


dm 

v d0 


(2.9) 


where and C pL are the density and specific heat of the 
drop liquid, T^ and T w are temperatures of the drop and free 
stream gas, A Is the heat of vaporization of the liquid, m 
is the droplet mass and 0 is the independent variable time. 
The vaporization rate, is given by Eq. (2.10). which follows 
from Eq, (2. 4) 


dm 

de 


" a 2 m a 


^Ao - 




(1 


( 2 . 10 ) 



where D, the droplet diameter, is. defined by Eq. (2.11) 



( 2 . 11 ) 


The mass transfer coefficients independent- of mass 
transfer rate are determined by the Ranz-Marshall correlation^ (46) 

Nu = 2 + 0.6 Re 1/2 Pr 1/3 (2.12) 


Nii^b = 2 + 0 . 6 Re 


Sc 1/3 


(2.13) 


where Reynolds, Prandtl and Schmidt numbers are designated 

by Re, Pr and Sc, respectively. The Nusselt, Nu, and 

i 

Sherwood, Nu A g, numbers are defined by Eqs. (2.14) and 
(2.15). 


Nu = 


hD 


(2.14) 


Nu 


k .DM 

A 


AB 


w 


(2.15) 


AB 


Average film properties of thermal conductivity (k), molecu- 
lar weight (M) , density (p) and binary diffusivity (P^g) 
are used. 

Combining Eqs. (2.4), (2.10) and (2.15), the result 
for drop mass is obtained. 



23 


Similarly, Eqs. (2.5), (2.7), (2.9) and (2.14) yield Eq. 
(2.17) for the drop temperature time derivative. 


dT 


L _ 


6 ' 


a* " Pt o T d 2 

pL 


k Nu z (T -T t ) . .. 

00 L ' A dm 

z n ~ ttD d0 


(2.17) 


Three auxiliary equations are required. 


_ _ P ^AB Nu AB M A ln(1 + R AB^ fr , H p v 
k ft Nu O + V P A R P B> 


(2.18) 


D 


■ &r 


n r = - 


p M A (1 


- y Ao } 


P L RT L 


(2.19) 


( 2 . 20 ) 


Eq. (2.20) is derived by assuming that the molar flux of 
species B is equal to the concentration of B at the drop 
interface times the velocity of liquid species A. 

Eqs. (2.16) and (2.17) constitute a coupled set of 
first order differential equations. The equations, with 
their initial conditions, were solved numerically by a 
fourth order Runge-Kutta method which has an automatic 
step size adjustment (56). It takes less than one minute 
to solve a three second drop history on the UNIVAC 1108 
computer. The film properties p, ^ AB , U, k, C , and M 
are determined at average film temperature and composition 



24 


values, and the Nusselt and Sherwood numbers are based 
upon average film values. The gas mixture component speci- 
fic heats, Cp A and C^g, are evaluated at the average film 
temperature. The properties p L and C ^ refer to the liquid 
at the liquid temperature. 

The steady state result is obtained by setting the 
derivative in Eq. (2.17) to zero. After substitution for 
z and rearrangement, 


m m _ * (e Z - 1) 

« ~ i L c . + N„ c _ 
pA R pB 


( 2 . 21 ) 


The equation for the steady state mass transfer rate is 
given by Eq. (2.16). The terms on the right hand side of 
the equation are functions of T^ and Tg, so it is necessary 
to solve by iteration. Since the right side is less sensi- 
tive to T , it was found most convenient to treat T as 
the dependent variable and T^ as the independent variable. 

A digital computer code for the IBM 1620 computer was 
written in FORGO language to perform steady state calcula- 
tions . 


Mixture Critical State Transition 

In the description of the physical process of droplet 
vaporization, the hypothesis of a transition from a vapori- 
zation process occuring at a constant plateau temperature 
to a steadily increasing one was advanced. The key to this 
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transition is the state of the gas mixture in the film 
surrounding the drop at the drop-film interface. If the 
gas mixture at this point exceeds the critical temperature 
corresponding to the total pressure in the gas mixture, an 
unsteady droplet temperature history results. Manrique 
and Borman's calculations (36) illustrate this phenomenon. 
Consequently, an estimate of the critical mixing line for 
the gas mixture is necessary to predict the low limit of 
this regime. 

A rigorous solution to this problem consists of 
solving the problem of thermodynamic equilibrium for the 
vapor-liquid interface. But in this case of the critical 
mixing line, one has to find the locus of the critical 
points of the y A (P, T) functions. Chueh and Prausnitz (13) 
in a recent article presented a method of estimating the 
critical mixing lines for non-polar binary mixtures. This 
method was adopted and programmed for computing the rela- 
tionships for the n-heptane-air and Fr-13-air systems. 

A brief outline of the method is given below. 

Chueh and Prausnitz proposed that the critical tempera' 
ture and critical volume be computed as a quadratic func- 
tion of surface fraction, 0^, 

a _ v v 2/3// „ 2/3 . v v 2/3\ 

®A “ y A V cA / ( y A V cA + y B cB J 


( 2 . 22 ) 
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where 

v ci - critical volume of component i, 

= mole fraction of component i. 

Parameters are introduced which are coefficients of the 
quadratic terms. They present correlations ' for these para- 
meters which are functions of the pure component critical 
properties. 

Because of the strongly nonlinear dependence of the 
critical pressure on composition, a different procedure is 
presented for computing the mixture critical pressure. 

The correlations for critical temperature and critical vol- 
ume are used in conjunction with a slightly modified ver- 
sion of the Redlich-Kwong equation of state. They compared 
computed critical pressures for 36 binary systems with data 
and found a mean average deviation of 3.6# The average 
deviations of experimental and fitted critical temperatures 
are found to be 0.4$ for 65 binary systems. Their scheme 
for determining the mixture critical volume was found to 
correlate the data for 25 systems within an average devia- 
tion of 1.9$. 


Property Non-Idealities 

The low pressure model departs from a realistic repre 
sentation of droplet vaporization as the system pressure 
increases. Considerations of the critical state of the 
gas mixture were presented in the previous section. This 
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section is concerned with the variations in the properties 
which enter into the film theory model. High pressure 
non-ideal effects are treated first and are followed by a 
discussion of the spatial variation of properties through 
the boundary layer. The implications of the use of pro- 
perties averaged over the film are assessed. 

High Pressure Effects 

The properties of the gas mixture, y Ao s X, k, C , P^, 

U, and p, all appear in the film theory results given in 
Eqs. (2.16) and (2.17). The first three, which are corrected 
in the high pressure application of the film theory, are 
discussed individually in subsequent parts of this section. 
The last four properties are used in their uncorrected form 
in the high pressure theory. Corrections for them were 
considered but discarded for various reasons. 

Calculations were made for a Fr-13 drop vaporizing at 
high pressure within the range of experimental conditions 
reported in Chapter V. Then arbitrary variations in the 
various mixture properties were introduced, and the sensi- 
tivity of the calculated plateau temperature and instantan- 
eous vaporization rate, to these changes were assessed. 

Negligible effects due to 10 % variations in mixture 

specific heat, C , viscosity, u, and density, p, were 

1 ? 

found. Therefore, high pressure corrections were not intro- 
duced for these properties. Since these values enter into 
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the dimensionless numbers which appear fractionally expon- 
entiated in the heat and mass transfer correlations , the 
insensitivity is not too surprising. 

However, a 10# variation in the binary diffusion coeffi- 
cient, produced about a 10# change in the steady state 

drop-air temperature difference and a 2 % change in vapori- 
zation rate. No data exists for binary diffusion coeffi- 
cients in Fr-13-air mixtures. The diffusivity values used 
were calculated by means of a corresponding state correla- 
tion which is estimated to be accurate only within approxi- 
mately 10# (25). Consequently, high pressure corrections 
to the binary diffusion coefficients (5*0 were omitted. 


Droplet-Film Interface Mole Fraction 

The determination of the mole fraction, y AQ , at the 
droplet-film interface is a vapor-liquid equilibrium pro- 
blem. In the case of a pure liquid drop vaporizing into a 
pure inert gas environment, there are two components and 
two phases. Designating the mole fraction in the liquid 
and vapor phases by x Aq and y AQ , respectively, of the four 
variable factors (y A , x A q , P, T) there are two degrees of 
freedom for the system according to Gibbs Phase Rule ( 4 0 ) . 

A common and simple approach is to use the saturation 
pressure relationship for the vaporizing component. 


^Ao ~ ^sat^total 


(2.23) 
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At total pressures low compared to the critical pressure 
of A, this usually gives a satisfactory approximation. 

This approximation is utilized in the low pressure film 
theory . 

At thermodynamic equilibrium the chemical potential 
of any given component in all phases of a multicomponent 
system is identical (40). This basic approach or its equi- 
valent of determining the mole fractions which give equi- 
valent fugacities in the liquid and vapor phases for each 
component was used by Manrique and Borman (36) in their 
analysis of the C.Og—Ng system. They specified the tempera- 
ture and total pressure of the system and determined the 
fugacity of the components as a function of component mole 
fractions using the Redlieh-Kwong equation of state for the 
mixture. In this manner they obtained a map of equilibrium 
mole fraction, y^ and x^, for the range of atmospheric to 
the gas mixture critical pressure for a given temperature. 

Manrique and Borman demonstrated their approach by 
making a satisfactory comparison with data for the C0 2 -N 2 
system. However, at the same time their approach is time 
consuming, particularly the algebra involved in obtaining 
simultaneous solutions to the component fugacities in both 
phases and satisfying species continuity in each phase. 

Hougen, Watson and Ragatz (23) describe another 
approach. It is applicable to the region near the critical 
point, and was used for the n-heptane-air calculations 



30 


given in Chapter V. This convergence pressure technique 
has been developed for hydrocarbon mixtures and extensive 
charts (39) have been prepared which give the vapor liquid 
equilibrium ratio, K, which is a function of P , T r , and 
the convergence pressure, P conv for each hydrocarbon, nitro- 
gen and carbon dioxide. The convergence pressure corres- 
ponds to the mixture critical pressure. Hougen, Watson 
and Ragatz indicated that this semi-empirical method re- 
quires judgment, but it yields generally acceptable results 
at pressures less than 80$ of the convergence pressure. 

A fourth method which was considered for use with the 
F r _13_ai r calculations and finally selected upon the basis 
of tests with data is summarized below. The method is a 
variation of the one outlined by Reid and Sherwood (47) 
for calculating phase-equilibrium constants, K^, 

K i = y i /x i (2 - 24) 


where the subscript o denoting the drop-film interface is 
dropped for convenience. 


K i “ y iL (f il/ P)/ ^i 


(2.25) 


where 

f° is the fugacity. of pure i at T and P as a liquid. 

XJLf 

If P > P gat , the saturation pressure at T, f|? L is estimated 
from the fugacity. of i saturated at T using the Poynting 
correction. If P < P sat , the saturated value of f? L 
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corresponding to T is used. For cases where T > T , 

c 

usually the case for nitrogen and air*, the vapor pressure 
relationship is extrapolated to a. fictitious value of 

The vapor phase fugacity coefficient for component i 
in the gas mixture, <j>^, is determined from the Redlich- 
Kwong equation of state as shown by Chao and Seeder (10). 

The liquid mixture is treated 
and the activity coefficient, Yj^* 

Scatchard-Hildebrand equation 

in y 1L * V* (6 1 - 6) 2 /RT 

where 

7 ■ 2 x i v2 

i i 

and 

v ± ■ ( v u )i <5 ' 7 + 3 V 

The molar liquid volume, V^, and the solubility parameter, 
6^, have been treated as empirical constants by various 
authors. Cavett (2.9). reviewed the literature and selected 
and tabulated values, for many components. Cavett's values 
for Ng, Og , COg , and n-heptane were utilized in test calcu- 
lations which are summarized below. 

In the case of Fr-12 and Fr-13 which were not listed 


as a regular solution 
is calculated by the 

( 2 . 26 ) 

V* (2.27) 


( 2 . 28 ) 
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by Cavett, (25° C) was determined by extrapolating the 
equation given by Prausnitz et_ al. (43). 

V ± = a + bT (2.29) 


Two values of the molal liquid volume at temperatures re- 
moved from the critical were substituted in order to evalu 
ate the constant in Eq. (2.29). Then (V^)^ was derived 
from Eq. (2.28) by substituting the 25° G value of V^. 

Cavett* s method (8) of determining <S^ was applied to 

give 


6 


i 




RT 



(2.30) 


y * 

where (H - H) is the enthalpy derivation for the liquid 
at 25° C, V is the (25° C) previously found and P is 
the saturation value at 25° C. Table I contains the con- 
stants derived for Fr-12 and Pr-13. 


Table I 


Calculated Constants for Liquid Phase Activity Coefficients 

1/2 


Component i 


(V J 


l cm 3 ^ 6. 
o> / l \gm mole / 1 \cm7 


Fr-12 

11.53 

6.96 

Fr-13 

11.86 

4.19 

Air 

2.60 

2.97 


An iterative procedure was developed to calculate the 
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equilibrium mole fractions for a two component mixture at 
a fixed temperature and pressure. Initially, K A was assumed 
to be equal to P ga t /P, Kg = °°, and x A = 1. Then, 

X B = (1 ~ X A K A )/K B (2.3D 

X A = 1 “ X B (2.32) 


y 


A 


"A 


K, 


(2.33) 


y B = 1 " y A (2.3 1 *) 

These four equations were solved and K A and K g were calcu- 
lated from Eq. (2.25). These calculated values were used 
as a second estimate and the procedure repeated until the 
fractional change in K A between successive iteration was 
less than 0.005. 

In order to avoid the difficulties of treating a three 
component mixture, air was assumed to be a single component. 
The constants in the liquid phase activity coefficient cal- 
culation were determined by mole fraction weighting of Ng 
and 0 2 values for the air composition. They are given in 
Table I. 

The results of a comparison of calculated and measured 
mole fractions for the COg-Ng system are shown on Fig. 1 
(60). The calculated gas phase values agree within 15 % at 
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P r % 1.5 and T r = 0.9 with respect to CC^. 

Three methods of calculating equilibrium mole frac- 
tions for n-heptane in air are compared with data ( 11 ) on 
Pig. 2. A constant pressure of P ^ 3 and a subcritical 
range of temperatures is covered. The calculated phase equi- 
librium constant method gives the best agreement. 

Pr-12-air system calculations are compared with the 
data of Parmelee (4l) in Pig. 3. The agreement over the 
range of data, up to P r of 0 . 3 , is good. 

The Chemical Abstracts from 1955-1969 were searched 
for vapor-liquid equilibrium data of any sort for Fr-13. 

The compilation of Chu et al. (12) published in 1956 was 
consulted in addition. No data was found. This negative 
result, of course, provided the motivation for developing 
and testing a method of computing the equilibrium mole fractions 
of binary mixtures at high pressures. 

Only the case of thermodynamic equilibrium at the 
interface has been considered to this point. It seems 
possible that the concentration of the continuous phase 
gas in the liquid may not reach the equilibrium value par- 
ticularly if the vaporization rates are rapid. A limiting 
value would be the case of no absorption in the vaporizing 
liquid. The phase equilibrium constants and, thence, the 
interface mole fractions for the limiting case of no absorp- 
tion can be easily calculated by the method outlined if the 
liquid phase activity coefficient, in Eq. (2.25) and 



N-Heptane Mole Fraction in Gas Phase 



g. 2. Equilibrium' Gas Phase N-Heptane Mole Fraction in 
N-Heptane-Nltrogen Two Phase Mixture at 100 Atm 
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in Eq. (2.32) are unity. The results of calculations 
employing this assumption for Pr-13 drops vaporizing in 
air are discussed in Chapter V. 

Latent Heat of Vaporization 

Consider a binary mixture in vapor-liquid equilibrium 
such as is normally assumed to exist at a vaporizing drop- 
film interface. The latent heat of vaporization of A, 
is defined as the difference between the partial molar 
enthalpy of A in the vapor and the enthalpy of A in the 
liquid. 

At low total pressures the assumption of substituting 
the heat of vaporization of pure A at T * T . is commonly 
used. This assumption is made in the low pressure film 
theory calculations. However, this approach will lead to 
substantial errors at total pressures in the vicinity of or 
greater than critical pressure of component A. A method 
of accounting for these high pressure effects has been de- 
veloped by Manrique and Borman (36). Their method which is 
applied to the determination of partial molal enthalpies 
of gas mixture components for use in the high pressure film 
theory calculations is outlined below. 

The Redlich-Kwong equation of state is applied to a 
gaseous mixture. The equation for the mixture compressi- 
bility factor is rearranged into the. form of a cubic and 
solved by the Newton-Raphson method. The enthalpy deviation 
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from the ideal state is computed from the temperature par- 
tial derivative of the log fugacity partial pressure ratio. 

Reid and Sherwood (47) have tabulated generalized 
enthalpy deviation as a function of reduced pressure, tem- 
perature and critical compressibility factor. Their gener- 
alized results based upon an integrated three parameter 
equation of state were extended to liquid phase enthalpy 
deviations by the use of a generalized correlation of latent 
heat of vaporization. The Reid and Sherwood values were 
used to find the liquid phase enthalpy deviations. 

Thus, the latent heat of vaporization in the high 
pressure region was computed as the difference of gas and 
liquid phase enthalpy deviations. Typical computed values 
for the Fr-13-air system are compared with the pure Fr-13 
values in Fig. 4. 

Thermal Conductivity 

Lenoir, Junk and Comings (34) developed a corresponding- 

states-type correlation for gases at high pressure. Their 

correlation is based upon nitrogen, methane, argon, ethane, 

ethylene, and carbon dioxide data over a large pressure range. 

§ 

Thermal conductivity, ratio, k , is correlated as a function 
of reduced pressure and temperature. It is defined as the 
ratio of conductivity at any pressure divided by the ther- 
mal conductivity at atmospheric pressure. 

Values of thermal conductivity for the pure components 




of the film gas mixture, e.g. , Fr-13 and air, are computed 
at the average film temperature. Atmospheric values are 
utilized. The low pressure values, k^ , are combined by 
means of the Mason-Saxena method which is given by Bird, 
Stewart and Lightfoot (4). They are weighted by mean film 
values of component mole fractions. 


where 
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(2.37) 


k^ is the mixture thermal conductivity at atmospheric 
pressure, k is the corrected value of mixture thermal con- 
ductivity, and y^ are the molecular weights and the vis- 
cosity of species i, and y^ is the film average mole frac- 
tion of species i. 

In order to correct this value for high pressure, the 
correlation of Lenoir, Junk and Comings is utilized with the 
pseudocritical properties determined by mole fraction aver- 
ages of the pure component critical temperature and pressure 
values. Then the reduced pressure and temperature are 



based on the total pressure and average film temperature 
values. 

Operationally, the values from the Lenoir, Junk and 
Comings correlation at various reduced temperatures at a 
specific reduced pressure were fitted by a third order 
least squares routine. This derived function was then 
utilized for Fr-13 drop vaporization calculations at a 
particular pressure. The reduced pressure could be fixed 
because fortunately the critical pressures of air and Fr-13 
are almost identical, and therefore, the pseudocritical 
mixture pressure is independent of mixture composition. 

Property Variation Through the Boundary Layer 

A major uncertainty in the application of film theory 
to droplet vaporization problems at high pressures and 
temperatures is the variation of properties through the 
boundary layer. Consideration was given to including a 
correction term or providing an averaging technique to 
account for this. This approach was abandoned after tech- 
niques were investigated. The techniques considered were 
found to be complicated and uncertain since no comparisons 
with data were given. Nevertheless, the techniques sur- 
veyed are interesting, and they are briefly summarized. 

The matched asymptotic series (1), reference state (32), 
successive iteration (52), and stagnation point methods 


are discussed. 
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Acrivos (1) showed that by piecing together exact 
asymptotic solutions for the cases of Prandtl or Schmidt 
numbers equal to zero, Nu Q , and infinity, Nu^, a simple 
expression for the Nusselt number in laminar boundary layer 
flow can be obtained. In the case of the constant property 
solution Acrivos demonstrated that the Nusselt number 
divided by the square root of the Reynolds number reduced 
by the asymptotic result for Prandtl number equal to zero 
is a simple function of NUq/Nu^. Comparison of exact numer- 
ical solutions indicated that his results were almost insen- 
sitive to geometry. 

He generalized this result to the case of variable 
properties. However, in order to solve for the Nusselt 
number, it is necessary to solve two separate non-linear 
integro-differential equations by expansion and summing of 
an asymptotic series. Derivatives and integrals of the 
properties have to be evaluated, and the algebra is very 
complicated. His approach is not easily adaptable to a 
digital computer. It would be necessary to further gener- 
alize his result for the case of simultaneous heat and mass 
transfer which would add further complication. Acrivos 
admits that there may be some question about the independence 
of the variable property result from geometry. 

A method of reference states was proposed by Knuth (32). 
He developed reference state prescriptions, by. a combination 
of exact calculations for laminar couette flow and the 
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analysis of few calculated results, of boundary layer prob- 
lems. There is substantial uncertainty in applying his 
results to correlations for over all transfer numbers to 
geometries other than flat plate particularly where separa- 
tion occurs. Furthermore, he claims that the composition 
reference state reduces to a linear average in all cases 
except where the ratio of molecular weights of the two 
species exceeds a factor of ten. The deviation of his 
reference states from those obtained by linear averages of 
temperature and composition were very small for a typical 
Fr-13 case. 

In the case of laminar boundary layer flow over a 
flat plate Schuh (52) performed calculations of variable 
properties. His method is quite simple and appears to be 
well suited for a digital computer. He starts with the 
constant property Blasius and Pohlhausen solution for the 
velocity and temperature profiles, determines the velocity 
profile with the viscosity and density variation from the 
initial temperature variation, and finally the new tempera- 
ture profile from the new. velocity and initial density 
profiles. It is claimed that this procedure converges after 
four repetitions. 

Schuh’ s calculations were applied to air with variable 
properties. He found that although the boundary layers 
thickened, there was no appreciable change in shearing 
stress or heat transfer coefficient over the constant 



property case. The uncertainty of applying this method to 
correcting heat and mass transfer coefficients for spheres 
exists since the flat plate equation is used. 

The idea of applying Schuh* s • successive iteration 
approach to the forward stagnation point on a two or three 
dimensional body of revolution naturally arises. However, 
there is again the question of applicability of a correc- 
tion for a single point to an average value determined 
empirically. Intuitively, one would expect better results 
than in applying the flat plate or plane two dimensional 
result directly. 



III. EXPERIMENTAL APPARATUS 


The experimental apparatus provides a means of intro- 
ducing and supporting a droplet in a stream of heated high 
pressure air and obtaining a size and temperature history 
of the evaporating droplet. J. A. Rieart-Lowe (48) and 
L. D. Alexander (2) have previously described parts of 
this apparatus. A complete description of the rig which 
is shown in a photograph. Pig. 5, is given here. 

Rig Components 

The rig is conveniently subdivided into a heater, test 
section and probe assembly. These components are described 
in order. 

A cross sectional view of the heater is shown in Fig. 
6. It consists of a cylindrical flanged pressure vessel 
which contains a packed bed heated by four semi-circular 
pipe heaters. The heater is "wye" connected with a 220 v 
three phase a-c source which provides a maximum total power 
of 4.6 kw. 

The air flow passes up from the inlet nozzle over the 
top of the radiation shield, down between the radiation 
shield and heater units, and up through the heater bed to 
the test section. The structural and thermal design of the 
heater are discussed by Ricart-Lowe (48). 

The test section consists of a 12 in. high rectangular 
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Pig. 5. High Pressure Suspended Droplet Vaporization Apparatus 
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block of stainless steel with a 6 in. square cross section. 
Pig. 7 shows a cross sectional view. Two window assemblies, 
a probe assembly which is described below, an inlet calming 
section and nozzle, numerous penetrations, and a top flange 
complete the list of components of the test section. 

The calming section, which is composed of a series of 
screens and spacers, tends to break up turbulence and flat- 
ten the air stream temperature profile. A 1 in. diameter 
ASME nozzle (3) was utilized. A teflon sleeve with an in- 
side diameter of 1 in. fits upon the top or throat of the 
nozzle in order to isolate the region around the probe tip 
from free convection currents. 

Penetrations in the test section are provided for a 
pressure gauge line (1/2 NPT) and the probe thermocouple 
leads (3/8 NPT). The probe assembly and two window sub- 
assemblies bolt to the sides of the test section. The test 
section is clamped by means of bolts between the heater 
upper flange and another blank flange which rests on top of . 
the test section. A detailed discussion of the test section 
is given by Ricart-Lowe (48). 

The probe which introduces and supports the droplet 
is an independent assembly which bolts to the side of the 
test section. It .consists of a housing (not shown) and a 
brazed subassembly, of tubing shown in Pig. 8. The sub- 
assembly which is joined to the housing juts, out from the 
housing flange into the test section where it makes a 
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Fig. 8. Cross-Section of Probe Subassembly 
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connection with the inlet cooling water line. 

The coolant enters the end of the 1/4 in. diameter 
tube of the probe subassembly where it is sealed by an 
O-ring to a fitting on the inside wall of the test section. 
The coolant flow through the tube around the probe tip 
and out the annular space between the 18 gauge test liquid 
supply line and the 1/4 in. tube. 

The test liquid flows from the bottom of the conden- 
sor-reservoir through the activating valve {# 1 on Pig. 9), 
through the 1/8 in. diameter stainless steel tubing to the 
probe housing where it passes through a pressure seal in 
the probe housing. The test liquid then leaves the probe 
housing and flows into an 18 gauge stainless steel hypo- 
dermic tubing part of which is shown on Fig. 8. 

The test fluid passes down around the teflon insulated 
thermocouple lead wires to the 0.029 in. diameter teflon 
probe tip. The section of 14 gauge tubing is sealed at 
the top with epoxy glue. The probe tip is also joined to 
the tube with epoxy. The thermocouple lead wires are 
sealed by epoxy to their spaghetti insulation about an 
inch above their outlet from the tubing. The test liquid 
is forced out the teflon probe tip where it can flow down 
the lead wires and form droplets on the thermocouple bead. 

A teflon collar was placed around the probe thermo- 
couple leads immediately above the junction in the Fr-13 
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tests. This was. found necessary in order to support the 
droplets for the tests. Alexander (2). gives a discussion 
of the probe details. 


Plow Arrangement 

The flow diagram. Pig. 9, shows the functional arrange- 
ment of the rig components. Paths of the vaporizing medium, 
water coolant and air are included in the diagram. 

The vaporizing medium is introduced into the Fr-22 
cooled condensor-reservoir through valve #2. It is fed 
into the probe by activating valve #1 and utilizing a differ- 
ence in pressure between the test section and reservoir. 

The regulating values, either high or low pressure, 
maintain the air pressure in the rig. Air acts as a pres- 
surizing agent for the probe coolant and the condensor- 
reservoir as well as a heat transfer medium for droplet 
vaporization in the test section. Air is available from 
either a 2000 psi reservoir (dried by cooling to 32° F at 
1800 psia) or a house air line at 75 psi. Air flow through 
the rig is regulated by choking the flow from the rig out- 
let at valve #6. 

Water cools the', probe supply line. The coolant is 
maintained at the .system pressure by a wash bottle arrange- 
ment which has to be refilled about every 2 hours. Cooling 
coils in an ice bath cool the water as it flows from the 


bottle to the test section. 









Instrumentation 


Instrumentation is provided to. obtain droplet tempera- 
ture and size histories, air flow rate, system pressure 
and various operating temperatures and pressures for per- 
formance maintenance, safety and control of the rig. The 
flow diagram, Pig. 9, shows the location of the temperature 
and pressure sensors. 

The probe temperature and the gas meter inlet air 
temperatures are recorded continuously on a brush recorder. 

The other temperature signals are recorded on a multipoint 
chart recorder. 

Either of two precision pressure gauges are used to 
determine system pressure. The ranges are 0-750 psia and 
0-1500 psia, respectively; the accuracy is 0.1# of full 
scale reading for each. 

The volumetric air flow through the test section is 

O 

determined by the use of a 1800 ft J /hr capacity gas meter 
and stop watch. 

The optical system used to obtain droplet size histories 
is shown schematically in Fig. 10. This arrangement pro- 
duces a shadowgraph of. the probe tip on a ground glass screen. 
Light from a point, source zirconium-arc lamp is focused about 
an inch in front of the thermocouple, and the shadow image 
is. collimated by a single lens upon the ground glass. The 
magnification usually, used is about six. Linear magnifica- 
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Fig. 10. Schematic Diagram' of Optical System 
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tion of objects at the probe plane was established by two 
methods. Optical bench mock-ups. of. the setup using still 
photographs of a reticle grid as well as actual objects 
placed in the probe position of the rig demonstrated line- 
arity. An image of a timer is projected on a mirror in 
the camera object plane. 

16 mm movie films are taken of the ground glass with 
a Cine-Kodak Special II movie camera. A 50 mm f/1.6 lens 
which focusses to 13 in. distance is used. 4 X reversible 
film with normal development procedure, lens wide open, 

13 in. object to film plane distance and 64 fps speed were 
found to give adequate size histories from which measure- 
ments could be taken. 

A switching circuit was built to synchronize the 
temperature and size histories. When a switch is thrown, 
a timing mark is produced on a channel of the temperature 
recorder, and the timer which appears on the movie record 
is started. 


Rig Operation 

A description of. the operating procedure for data 
gathering is given in this section. 

The heater is. turned on and house air is flowed through 
the rig until the entire assembly is warmed up. This takes 
about six hours. During the start up period the instrumen- 
tation is warmed up and calibrated, the coolant cylinder 
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and ice bath are filled, and the optical system is aligned. 

The test fluid is placed into the coils of the con- 
densor. After the rig has achieved steady state, the air 
fluid reservoir is pressurized to about 75 psi above the 
desired operating pressure. The probe coolant is started, 
and the system pressure and the air flow rate are adjusted 
to the test values. 

At this stage the test is initiated as follows. The 
test liquid can be made to drip on the probe thermocouple 
by manipulating valve #1. The temperature recorder is 
started. When a stable drop is residing on the probe thermo 
couple, the camera is started and the "start" switch of 
the synchronization circuit is activated. The camera is 
stopped when the drop disappears and the temperature recor- 
der remains on until the air temperature measured by the 
probe sensor recovers its initial value. Normally, the air 
flow is stopped to conserve air and a cursory analysis of 
the temperature history is made before proceeding to the 
next droplet history. 
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IV. EXPERIMENTAL UNCERTAINTIES. AND DATA REDUCTION 

In the opening section of this chapter experimental 
problems which introduce errors are discussed. The second 
section deals with the data reduction techniques which were 
employed. Uncertainties in the experimental results and 
their implications regarding the calculated results are 
treated in the final section. 

Sources of Experimental Error 

There are several problems which attend high pressure 
suspended droplet vaporization experiments. The combination 
of forced and free convection, reduced surface tension, and 
non-convection heat losses are among the problems; they 
are discussed below. 

In mixed or combined convection problems, the important 
parameter is Gr/Re , the ratio of the Grashof number to 
Reynolds number squared (30). In the situation where the 
characteristic length, forced flow velocity, fluid tempera- 

p 

ture and temperature difference remain constant, Gr/Re 

would be independent of pressure. If these conditions did 

indeed prevail for droplet vaporization, the importance of 

2 

combined convection, which is measured by Gr/Re , would not 
increase as the pressure is increased. 

Experimentally it was found that the air velocity in 
the test section had to be reduced as the pressure was 
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increased in order to maintain stable drops. Since the 

2 

other parameters in Gr/Re remained approximately constant, 

2 

Gr/Re increased with increasing pressure. A consideration 
of the critical Weber number criteria for drop stability 
explains this experimental fact. 

We = U - ^ - D - P - (4.1) 

As shown in Eq. (4.1) the square of the velocity, u, corres- 
ponding to the critical value of the Weber number is propor- 
tional to drop surface tension, a, and inversely proportional 
to the gas density, p. Through the combined effect of gas 
density which increases with pressure and drop surface ten- 
sion which decreases with pressure (47), the maximum velocity 
that results in a stable droplet decreases with increasing 
pressure. 

In the region of combined free and forced convection 

2 

heat transfer to spheres, Gr/Re > 0.3, the heat transfer 
rates for opposed flow are lower than that predicted by pure 
forced convection (59). However, for the case of simultaneous 
heat and mass transfer to spheres in opposed flow, the oppo- 
site effect has been observed although the region of combined 
convection, approximately 0.001 < Gr/Re 2 < 200 for = 1, 
is different (38). Although velocity components arise from 
both temperature and concentration differences, both are in 
the same direction, i.e., downward for cool, dense vapors. 

In the case of Fr-13 data, both are important but neither 
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predominates. This rather remarkable result is another 
illustration of the danger inherent in indiscriminate appli- 
cation of heat-mass transfer analogies. Since Narasimhan 
and Gauvin's experiments (38) were carried out at turbu- 
lence intensity levels less than \% > at Prandtl numbers of 
1, at steady state and over a limited range of variables, 
their results were only used interpretatively and not 
quantitatively. 

The upper limit on velocity at higher pressures also 

''' 

causes an increase in the free convection from the test 
section walls. This tends to disturb the uniform velocity 
distribution at the outlet of the nozzle. In the initial 
high pressure experiments fluctuations in air temperature 
at the probe tip led to two successive design changes in 
the apparatus. The first change was an enlargement of the 
flow nozzle from 1/2 to 1 in. diameter. The second, 
which finally eliminated the problem, was a lengthened 
nozzle. 

A second source of experimental error relates to 
reduced droplet surface tension. Since the droplet surface 
tension decreases with increasing pressure, droplet stabi- 
lity becomes a problem. This stability phenomenon not only 
limits the velocity which enhances free convection as dis- 
cussed above, but it also leads to troublesome droplet 
behavior, and it limits the droplet size obtainable. 

Glasstone (20) gives the following formula for the 
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apparent weight , W, of a drop which drips, from a fed tube 
of radius, r, 

W = 2irro(f) (4.2) 

where <j> is an empirically determined correction factor. 

In applying this criterion to estimate drop size, r would 
correspond to the probe tube tip radius. An increase in 
ambient pressure would decrease a which in turn would limit 
the drop weight or size obtainable. 

At the onset of instability, the droplet oscillates in 
a vertical direction alternately taking the classical ob- 
late and prolate shapes. When this problem occurs at high 
pressures, it either leads to loss of support (droplet falls) 
or damped droplet oscillations after sufficient size reduc- 
tion due to rapid vaporization. Occasionally droplet spin- 
ning occurred at high pressures. The effects of this be- 
havior on heat and mass transfer rates are difficult to 
assess, but it undoubtedly increased the average rates 
above what would be predicted by steady state correlations. 

Analysis of the movie film histories from an initial 
series of Fr-13 experiments showed that frequently the drop- 
let clung to the bottom half of the spherical thermocouple 
bead and did not envelop it as most liquids would at room 
temperature and atmospheric pressure. The exposure of the 
unwetted top portion of the bead provided an extraneous 
heat energy source to the droplet. The plateau temperatures 
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obtained in this series of tests were suspiciously high 
when compared to values predicted by film theory. This 
was later confirmed by Fr-13 data taken with a redesigned 
probe tip. The latter data are reported in Chapter V. 

There are two non- convection sources of heat transfer 
to a suspended drop which may contribute to experimental 
error. One is radiation from the test section walls, and 
the other is the thermocouple lead wires. An estimate of 
the fraction of the total heat rate transferred to the drop 
by radiation was made for a typical experimental condition. 

The radiative energy flux was considered to be composed 
of two components. The first from the test section walls 
which were assumed to be at 150° F (the outside wall is 
warm to the touch) was assumed to have emissivity and view 
factors of unity. A second component was due to the air 
calming section which can be seen through the nozzle from 
the drop location. It was assumed that the calming section 
was at 300° F, the apparatus maximum air temperature. A 
view factor was estimated by assuming that the geometry 
corresponded to two parallel disks (33), and that the effec- 
tive emissivity was unity. 

The computed radiative heat flux was found to be less 
than b% of the convective heat flux for an 0.033 in. diameter 
drop with a computed Nusselt number of 16.4 and air-drop 
temperature difference of 100° F. 

An estimate of the heat loss from the thermocouple lead 
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wires was also made. It was assumed that the 0.005 in. 
diameter lead wires acted as infinitely long one-dimensional 
fins. An average wire thermal conductivity of 10 BTU/hr- 
ft-°F was assumed which corresponds to the chromel-constantan 
lead wires. 

A temperature difference between the lead wire root 
and ambient gas of 180° F and a convective heat transfer 
coefficient of 230 BTU/hr-ft 2 -°F were used in the calcula- 
tion. The heat energy arriving from both lead wires com- 
puted in this manner was found to be less than 3 % of the 
convective heat energy found for the previously outlined 
case. Although the model assumed is crude, it is believed 
that the assumptions are conservative. Furthermore, for 
the Fr-13 data, 0.003 in. diameter lead wires were substi- 
tuted for the 0.005 in. wires assumed in this calculation. 

Data Reduction Techniques 

The important data from a drop vaporization run are 
droplet temperature history, size history, air temperature, 
air pressure and gas flow rate through the test section. 

In comparing the results of transient calculations with 
experimental droplet, histories, the calculated dependent 
variables are droplet temperature and size histories. Repre- 
sentative experimental histories for a Fr-13 drop run are 
illustrated below, and then the data reduction technique 


is described. 
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The droplet size histories were obtained on 16 mm 
movie film. Prints, of a Fr-13 film history are shown on 
Fig. 11. A photograph of the corresponding temperature 
history record for the Brush recorder is shown on Fig. 12. 

The conditions of the run were pressure = 47.8 atm (P = 
1.25), air temperature = 171.5° F, and air velocity 3.37 
in. /sec. The timer started about 32 frames before frame 0. 
The voltage step shown on the channel 1 of the temperature 
record coincides with the timer commencement. The marks at 
the top edge of the temperature record are second marks; 
the chart speed was 125 mm/sec. The droplet was being fed 
until frame 2. 

The temperature history is reduced by picking milli- 
volt values off the data record at frequent intervals and 
translating them to temperatures by means of a conversion 
chart. Each channel of the recorder was calibrated with a 
potentiometer over the entire range being used, and this 
procedure was repeated each time the range was changed. 

A rather elaborate procedure of droplet size data 
reduction was necessary because of the drop shape distortions 
The data shown in Fig. 11 illustrate the non-spherical drop 
silhouettes. After a particular frame of the film record 
was selected for size determination, it was projected on a 
large white piece of cardboard at a range of about 15 ft. 

The drop outline was traced on a piece of paper. Next, the 
paper was put on a glass-topped tracing table, and the drop 




Pig. 11. Photographic Drop History of Fr-13 Drop Vaporizing in Air at 47.8 atm abs 
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outline was traced on vellum paper. 

An axis of symmetry of the drop was found by folding 
the vellum sheet, so that it split the drop profile in 
half. This step required judgment since the drop half out- 
lines were not always coincident. For the drop profile 
shown on Fig. 13 about 12 # of the measured drop volume is 
contained in the region between the overlapping boundaries. 
8 % is positive and 4# is negative. One half of the drop 
profile was divided into a number of segments the same 
height as shown on Fig. 13* Then the drop volume and sur- 
face area were determined by assuming that the drop was 
composed of frustrum of a cone shaped segments with radii 
shown on Fig. 13- The drop volume, V, and surface area, 

A, were computed by Eqs. (4.3) and (4.4). 



n-1 is the number of increments, r are the radii, and h is 

n 

the segment height. The mean chord length, D ch , was cal- 
culated by Eq. (4.5) 


D oh = 6V/A 


(4.5) 



Drop outline 


1 



Pig. 13. Profile of Fr-13 Drop Vaporizing in Air at 
47.8 atm abs 
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The mean chord length is used in the calculations as a 
characteristic length in the heat and mass transfer corre- 
lations and to relate the droplet mass to area with Eq. 

(4.5). However, a further relation is required to relate 
D k to droplet volume. This relationship is discussed in 
the following paragraphs. 

Since the drops were not spherical, D Qh differed from 
the diameter of an equivalent volume sphere. The values 
computed from Eq. (4.5) were plotted versus volume, V, for 
a number of data points. This function which is shown on 
Fig. 14 is within a slight scatter, independent of tempera- 
ture and pressure. The third order least squares fit to 
this function and the diameter for an equivalent volume 
sphere are shown on the figure. The derived relationship 
between D ch and V was utilized in the droplet history calcu- 
lations presented in Chapter V. 

Droplet mass and temperature were the dependent variables 
in the transient calculations which are compared with measured 
droplet histories. At each time interval in the calculation, 
the droplet volume is determined from the mass, and then the 
D ch (V) derived relationship gives the diameter and surface 
area which are used to. calculate instantaneous time deriva- 
tives of droplet mass and temperature. In this fashion the 
experimental surface to volume relationship is used In ob- 
taining the calculated droplet histories. 
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Estimated Experimental Errors 

In the first section of this chapter sources of error 
which would be expected to introduce bias into the experi- 
mental results were discussed. Accurate predictions of 
these errors is very difficult, and about the only recourse 
available is to design the apparatus and experimental pro- 
cedure to minimize these effects. In addition to such errors 
there are, of course, experimental uncertainties. 

Experimental uncertainties include accuracy errors and 
precision errors. Accuracy errors are detected by calibra- 
tion procedures, the results of which are applied when the 
data is processed. This was done for the pressure, tempera- 
ture, gas flow rate and linear size measurements. Precision 
errors are due to random fluctuations of the instrument it- 
self in conjunction with the data reading process. An accu- 
rate estimate of precision errors requires an experimental 
determination of measurement error which, since it involves 
subsidiary experiments, is expensive and time consuming. 

In lieu of experimental determinations of precision 
errors, uncertainties can be readily estimated. Uncertainty 
is defined as "what we think the error would be if we could 
and did measure it by calibration” (51). An approximate 
value of uncertainty that is assumed is a plus or minus accu- 
racy range which includes 95 - 99 % of all readings of the 
instrument (51). If no accuracy range is given by the 
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instrument manufacturers, one-half the least count or scale 
reading is used. Thus, with these definitions and the major . 
assumption that the individual deviation populations are 
normally distributed, standard procedures for combining inde- 
pendent errors can be applied to obtain precision indices 
of the experimental results. Furthermore, the implications 
of these errors regarding the calculated results compared 
with experimental results in Chapter V can be estimated. 

Measurement Un cert aint ie s 

The following list of assumptions leads to a simple 
statistical technique ( 51 ) for estimating the uncertainty 
interval, R, in a result, R, due to a combination of errors, 
x, in variables, x. 

1. Estimates of uncertainty intervals are valid for 
the various variables. 

2. Measurement deviation populations are normally 
distributed. 

3. Errors are independent. 

The uncertainty interval for N variables is given by Eq. (4.6) 

x=l v / 

This technique is utilized to predict the uncertainty 
interval for the gas flow measurement, (/, the velocity in 
the test section, u, the temperatures of the drop, T^ , 
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air, T , and gas meter inlet, T , and the measured mass 

' oo J g ' 

transfer rate, The gas flow is determined by timing 

the interval for a given quantity of gas to flow through 
the meter as given by Eq. (4.7) 


V = | (4.7) 

"5 

where V is the total volume, typically 3 ft , and 6 is the 
stopwatch time, typically 1 min. Applying Eq. (4.6) 




(4.8) 


A measured temperature, T, consists of a thermocouple 
potential, v tc , multiplied by scale factors, K^, as given 
by Eq. (4.9) 

T " v tc K 1 K 2 K 3 (4 * 9 ? 

where the scale factors are the gain of the amplifier, the 
recorder scale factor, and the slope of the thermocouple 
calibration curve, respectively. Applying Eq. (4.6), the 
uncertainty interval for the temperature measurement is 
obtained as Eq. (4.10). 

(i) 2 =@ 2 4 2 /+(i§) 2 

The test section velocity, u, is computed from Eq. 


(4.11) 
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u 


C 1 /. T 

00 

T“t 
. g 


(4.11) 


where C is a constant which depends on the barometric pres- 
sure and test section nozzle size which are assumed to be 
accurately known, and T^ is the gas meter inlet temperature. 
Similar manipulations to Eq. (4.11) give the uncertainty 
interval for u. 



(4.12) 


The slope of the drop mass versus time gives the drop 
mass transfer rate, dm/d0. The liquid density, of course, 
has to be used to convert from drop volume to mass. In 
order to obtain an estimate of measured mass transfer rate 
uncertainty interval, dm/d0, it is assumed that the measured 
drop mass transfer rate is proportional to the cube of the 
drop diameter, D, and inversely proportional to drop life- 
time, A0. 



C D 3 
A0 


(4.13) 


where C is a constant which includes the liquid density. 
Operating on Eq. (4.13) with Eq. (4.6), Eq. (4.14) results. 



(4.14) 
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The uncertainty interval in droplet diameter, D, is 
treated in a subsequent paragraph and its. calculation is 
given in Eqs. (4.15) and (4.16). 

A summary of the uncertainty intervals for the funda- 
mental components of the experimental measured quantities 
along with their sources are contained in Table II. Table 
III contains a summary of the calculated uncertainty inter- 
vals for a typical Fr-13-air experimental test run. The 
typical values selected were P = 38.2 atm (P = 1), u = 
u = 3 in. /sec, T * 220° P, T * 70° F, D = 0.04 in., and 

s 

T l = -13° F. 

Table II 

Summary of Measurement Component Uncertainty Intervals 
Uncertainty, x 

Quantity, x (absolute) Source 

V 0.05 ft ^ 1/2 of minimum scale 

division 

0 0.005 min 1/2 of minimum scale 

division 

1 mv/v Manufacturer's specifications 

K„ 0.025 div/full 1/2 of minimum scale 

d range division 

K_ 1.5° R/v Slope of t.c. cal. curve 

3 for probe 

K_ 2.5° R/v Slope of t.c, cal. curve 

3 for g.m. temp. 

A0 0.01 sec 1/2 of minimum readable 

timer scale division 
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Table III 

Summary of Uncertainty Intervals Calculated for a Fr-13 Test 
Uncertainty Interval 


Quantity 

absolute 

% 

Comments 




P 

0 . 1 atm 

0.3 

Gauge mfg. 

spec . 


V 

0.05 ft 3 /min 

1.7 

Calculated by 

Eq. 

(4.8) 

u 

0.2 in. /sec 

4.4 

Calculated 

by 

Eq. 

(4.12) 

t l 

1.5° R 

3 

Calculated 

by 

Eq. 

(4.10) 

T 

g 

1° R 

2.7 

Calculated 

by 

Eq. 

(4.10) 

T 

OO 

5.6° R 

3 

Calculated 

by 

Eq. 

(4.10) 

D 

0.00082 in. 

2 

Calculated 

by 

Eq. 

(4.16) 

dm/d0 

1.3xl0 - ^ lb /sec 
m 

6.5 

Calculated 

by 

Eq. 

(4.14) 

The 

determination of the 

mean 

chord length < 

of 

a drop- 


let from film data was described in the previous section. 
The experimentally derived relationship between the droplet 
mean chord length and volume was shown on Fig. 14. The 
function was smoothed and fitted by a third order least 
squares fit through the 20 data points. The scatter in 
these data about the fitted function gives the deviation of 
the data from a mean value if the deviations are approxi- 
mately normally distributed. 

A test of the normality of the data is shown on Fig. 

15 where the data is plotted on probability coordinates. 
Since the data falls in approximately a straight line which 
passes through the center, the distribution is symmetrical 



Percentage of deviation at or belo’ 



Pig. 15. Distribution of Measured Mean Chord Lengths Prom 
Least Square Values, Fr-13 Drop Data 
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and approximately normal (51). As a result Eq. (4.15) can 
be utilized to compute the standard deviation of the dia- 
meter measurement. a D . 

N 

S U L " X)2 
1 = 1 ' 

N - 1 

0 ^ was calculated from the data and the fit is shown 
in Pig. 14 and found to be 5.76 x 10 - ^ in. Eq. (4.16) gives 
the relationship between the uncertainty interval, x, and 
a D for the normal distribution. 





(4.16) 


where P(x) and P(a) are the probabilities that a given devia- 
tion will fall within the uncertainty interval, + x, and 
+ a, respectively. The results for the diameter uncertainty 
are also given in Table III. 


Reflected Uncertainties in Calculational Comparisons 


Comparisons between droplet history data and calculated 
droplet histories are presented in the following chapter. 

The calculated histories require values of independent 
variables which correspond to the experimental conditions 
of the test. The major dependent values which are calcu- 
lated are the air-droplet temperature difference and the 
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droplet mass transfer rate. Uncertainties in the values 
of these independent variables, u, D, T^, and P, lead to 
uncertainties in the major dependent variables, T^ - T L 
and dm/de. It is possible to estimate uncertainty inter- 
vals in these calculated variables. 

Eq. (4.6) was stated in the previous section. It can 
also be applied to the steady state droplet calculations 
developed in Chapter II. The partial derivatives in Eq. 

(4.6) in this case are non-linear and complicated with 
respect to their dependence on x. Consequently, the 3R/3x 
values for R = (T - T t ) and dm/d0 were determined by per- 
forming computer calculations of the change in the dependent 
variable, R, for small changes in the independent variables, 
x, about a base value. 

The base conditions which are typical of the Fr-13-air 
data are: P = 38.2 atm (P^ = 1), u = 3 in. /sec, T ro * 220° F, 

D = 0.04 in., and T L = -13° F. The uncertainty intervals, 

+ x, for u, D, T, and P were given in Table III. A summary 
of the partial derivatives and the components of the computed 
dependent variable uncertainty intervals is given in Table 
IV. 

The final total values of the uncertainty intervals are: 

T =TTT = 4.85° R 
00 L 


= 2.1% of base case, values and 



Table IV 
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Stt = 1.02 X 10" 7 lb /sec 
do m 

= 5.2$ of base case value. 

What these results mean Is that for a typical case uncer- 
tainties in the measured values of the independent variables 

introduce uncertainties of + 5° R and + 10~ 7 lb /sec, res- 

- - m 

pectively, into the calculated temperature difference and 
mass transfer rates. These values can be translated to 
probable error which is the conventional method of reporting 
results. Probable error here is defined as the range in 
which one half of the reported values would be expected to 
lie. The probable error interval would be approximately 
50 # of the uncertainty interval, x, and the computed results 

O 

would be reported as T m - T^ + 2.5 R and dm/d9 + 5 x 10 
lb m /sec. 

The other significant result of this analysis is that 
a good test of theoretical model would be to see if it con- 
sistently predicted results which agreed within the sum of 
the experimental and calculated uncertainty intervals. This 
would mean drop temperatures which agreed to within 13° R 
and mass transfer rates to within 23$. 
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V. EXPERIMENTAL RESULTS AND COMPARISONS WITH THEORY 

This chapter contains a presentation of drop vapori- 
zation data for water, n-heptane and Freon-13 (Fr-13). The 
data are compared with film theory calculated results for 
plateau temperatures and vaporization rates. The Fr-13 
temperature histories which exhibited non-classical behavior 
are correlated with the critical point of the gas mixture 
at the drop interface. 

Water Drops Vaporizing in Air 

A few temperature data points were obtained for water 
drops vaporizing in air at pressures of 20,40 and 60 atm 
abs. The water drop experimentation was undertaken for the 
purpose of testing the rig performance, and the data were 
previously reported by Alexander (2). 

An initial series of tests was performed utilizing a 
copper-constantan probe thermocouple which had a 0.0185 in. 
diameter spherical bead and 0,005 in. diameter lead wires. 
For this series of tests the probe tip consisted of 22 
gauge stainless steel hypodermic tubing instead of the tef- 
lon sleeve (described in Chapter III) which was installed 
later. 

The copper-constantan probe thermocouple was replaced 
with an iron- const ant an one of the same dimensions and the 
stainless probe tip was replaced with the teflon tube. 
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After these alterations were made, a second series of tests 
were made. The results, of both tests are shown on Pigs. 

16, 17 and 18 for the test pressures of 20, 40 and 60 atm 
abs, respectively. The results of low pressure film theory 
calculations of plateau temperatures are also shown in the 
figures for comparison purposes. 

Because of the relatively high surface tension of 
water, the drops clung to the probe tip instead of the 
thermocouple junction as was the case for n-heptane and 
Fr-13 drops. Initially, the drops were approximately 0.1 
in. diameter, and they completely enveloped the thermocouple. 
Although the drop lifetimes were long, on the order of 
minutes, no mass transfer data were obtained. The plateau 
temperatures which were measured with a potentiometer were 
steady. No air flow conditions were recorded. 

The 20 and 40 atm pressure results show that the probe 
alterations reduced the indicated wet bulb temperatures 
by about 15° P, and the corresponding reduction was 20° F 
for the 60 atm data. Since the lead wires were immersed in 
the drop for both series of tests, it is surmised that the 
substitution of the teflon for the stainless probe tip 
reduced the drop heat gain by conduction sufficiently to 
account for the difference. 

The better data, which were obtained with the iron- 
constantan thermocouple arrangement, agree with the calcu- 
lated plateau temperatures to within 10° P for all points. 








Cu-Con t.c. 
Fe-Con t.c. 
Calculations 



Fig. 18. Plateau vs Air Temperatures for Water Drops Vaporizing in Air at 
60 atm abs 
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However, the calculated values are consistently lower. 

The agreement between the data and calculations seems to 
improve as the test pressure is increased although the 
data are too meager for a conclusive statement about the 
application of low pressure film theory to the data. The 
generally good agreement between the calculated and measured 
drop temperature was considered to be adequate for attesting 
satisfactory rig performance at high pressures. 

N-Heptane Drops Vaporizing in Air 

About 90 histories were recorded for n-heptane drops 
vaporizing in air at pressures of 1.5, 5, 10, 50 and 100 
atm abs. The data covered a range of reduced pressures of 
0.056 to 3.71 and reduced drop temperatures of 0.57 to 
O .69 with respect to the critical properties of pure n- 
heptane. A chromal-constantan probe thermocouple with a 
bead diameter of 0.018 in. diameter and 0.005 in. diameter 
lead wires was used to obtain all n-heptane temperature 
data. All of the drop temperature histories displayed the 
typical behavior where the temperature initially sharply 
approaches the plateau or wet bulb temperature, levels out 
at the plateau value and slightly before the thermocouple 
is dry, begins to rise exponentially to the air temperature. 

The minimum test section pressure which would produce 
adequate air velocities was about 1.5 atm abs. Consequently, 
it was not possible to make an exact comparison with 
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published vaporization data for n-heptane drops in air at 
atmospheric pressure. However, calculated plateau tempera- 
tures were compared with both calculated and measured 
values reported by Priem (44). The results of this com- 
parison of atmospheric data and calculations are shown on 
Pig. 19. The agreement between the two calculated values 
is excellent both in magnitude and trend. The values 
calculated by the steady state technique outlined in Chapter 
II are consistently about 5° F higher than the experimental 
points of Priem. This agreement is within the experimental 
accuracy reported by Priem for 1 atm data. 

A comparison between calculated plateau temperatures 
and measured values at the lowest system pressure investi- 
gated, 1.5 atm abs, is shown on Pig. 20. An air tempera- 
ture range of 110-250° P is covered by the data. The agree- 
ment between calculated and measured values is within 10° P 
with a maximum scatter about a mean experimental fit of 
+ 5° F. Calculated and measured plateau temperatures agree 
better at the higher air temperatures. A few data points 
taken on a different day from the bulk of the data shown 
on Pig. 20 are indicated by the triangles. The fact that 
these points fall generally in the range of the other data 
demonstrates self consistency of the data and experimental 
technique. 

The difference between plateau temperatures calculated 




Pig. 19. Comparison Between Calculated Plateau Temperatures 
and Data of Priem, n-heptane Drops in Air at 
1 atm abs 



100 150 200 250 300 
Air Temperature, °F 


Fig. 20. Plateau vs Air Temperatures for N-Heptane Drops 
Vaporizing In Air at 1.5 atm abs 
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by low pressure film theory and film theory with the inter- 
face mole fraction and heat of vaporization corrected for 
pressure is almost indistinguishable at 1.5 atm. The high 
pressure film theory utilized in- calculating n-heptane drop 
steady state temperatures and mass transfer rates included 
no corrections for surface regression or gas film thermal 
conductivity as in the case of the Fr-13 analysis. 

Fig. 21 shows a comparison between calculated and 
measured plateau temperatures for n-heptane drops vaporizing 
in air at 5 atm abs. In the range of data shown, both the 
high and low pressure models agree with the data to within 
5° F except at air temperatures above 250° F where both 
models predict plateau temperatures 8-10° F higher. The 
film theory calculations corrected for high pressure are 
2-3° F closer to the measured values than the uncorrected 
theory for air temperatures above 250° F. 

10 atm abs data is compared with high and low pressure 
film theory results on Fig. 22. The high and low pressure 
film theory results agree within 5° F over the range of the 
comparison. The lower plateau temperatures predicted by the 
high pressure theory fit the data better than the low pres- 
sure theory, although both theoretical results are in good 
agreement with the data. 

A comparison of calculated and measured plateau tempera- 
tures of n-heptane drops vaporizing in air at 50 atm abs 
is shown on Fig. 23. This pressure is almost twice the 




100 150 200 250 300 
Air Temperature, °F 


Fig. 21. Plateau vs Air Temperatures for N-Heptane Drops 
Vaporizing in Air at 5 atm abs 




Fig. 22. Plateau vs Air Temperatures for N-Heptane Drops 
Vaporizing in Air at 10 atm abs 
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critical pressure of n-heptane. The shape of the droplet 
temperature histories at this supercritical pressure was 
similar to the data at the lower air pressures. The results 
calculated by either theory agree within 10° P of the 
measured values. The low pressure film theory calculated 
values show better agreement with the data. 

Initially, all of the high pressure film theory calcu- 
lations were based upon n-heptane interface mole fraction 
estimates obtained by the convergence pressure technique 
discussed in Chapter II. When the mass transfer rate com- 
parisons were made, the high pressure film theory results 
consistently predicted mass transfer rates about 30% higher 
than the measured values for the 50 and 100 atm cases. 

This anomaly was subsequently traced to the high predicted 
n-heptane gas phase mole fraction at the drop-film inter- 
face. A comparison of the results of the convergence pres- 
sure calculation with measured values was shown on Fig. 2 
for n-heptane-nitrogen mixtures at 100 atm. The predicted 
value of the n-heptane gas phase mole fraction was about 
40% higher than the data for the highest temperature point 
(388° K = 240° F) . 

The problem of inaccurate calculated mole fractions 
was resolved by using the n-heptane-nitrogen vapor-liquid 
equilibrium data in lieu of calculated values for the film 
theory vaporization calculations at 50 and 100 atm. The 
mole fractions calculated by the convergence pressure 
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method at 100 atm were appoxlmately 20# higher than those 
calculated at 50 atm and the values were equal at 260° F 
for 100 and 50 atm. Consequently, the 100 atm vapor-liquid 
equilibrium data were used unaltered in the 50 atm calcu- 
lations. That this procedure is more satisfactory than 
using the calculated values for all pressures is borne out 
by the resulting mass transfer comparisons which are pre- 
sented in the latter part of the section. The use of 
accurate values for the interface mole fraction at 50 atm, 
if they were available, would be expected to produce higher 
calculated plateau temperatures than those shown on Fig. 23 
and hence, probably better agreement. 

Calculated and measured plateau temperature results 
at 100 atm abs are compared on Fig. 24. The temperatures 
predicted by steady state film theory corrected for high 
pressure agree within 5° F for all data points except the 
highest which is about 10° F high. The low pressure theory 
results agree well with the data too, although they appear 
to be consistently a few degrees higher than the average 
data values. 

Fig. 25 shows a summary of the plateau temperature 
data at all pressures. The lines representing the tempera- 
ture data are least squares fits of the data points. The 
slopes of the plateau versus air temperature functions 
steepen as the pressure increases, and they approach the 
plateau equal to air temperature line as an asymptote as 




Fig. 24. Plateau vs Air Temperature for N-Heptane Drops 
Vaporizing in Air at 100 atm abs 
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the pressure Increases. The knee in the plateau versus 
air temperature curves shifts in the direction of higher 
air temperature as the system pressure is increased. This 
effect of pressure for a given fuel is analogous to the 
effect of going from less to more volatile fuels at a given 
pressure (44). 

The data represented on Fig. 25 are shown cross plotted 
versus pressure on Fig. 2 6 for air temperatures of 100, 

175 and 250° F, respectively. The slope increases with 
increasing air temperature although the lines appear to 
flatten at the highest pressures. They would most likely 
approach the air temperature as a limit at very high pres- 
sures . 

The effect of pressure upon the two film theory predic- 
tions is represented on Fig. 27 where the predicted drop 
temperature minus the measured plateau temperature is plotted 
versus pressure for an air temperature of 175° F. The 
measured plateau temperatures are least square fitted values. 
Temperature errors are within + 6° F for both theories at 
all pressures. The lov[ pressure film theory error exhibits 
a tendency to increase with increasing pressure. 

Several individual droplet histories were analyzed for 
mass transfer rates in the manner described in Chapter IV. 

In order to obtain instantaneous mass transfer rates, the 
drop volumes at several different times during the history 
were plotted versus time. A tangent to the volume time 



220 r— Air Temperature 



Fig. 26. Summary of N-Heptane Temperature Data, Shown vs Pressure for Air 
Temperatures of 100, 175 and 250°F 



Priem's Data (44) 



Pig. 27. Effect of Pressure on Calculated Plateau Temperature Error for 
N-Heptane Drops Vaporizing in Air at 175°F 



103 


curve was measured at the particular point desired. Care 
was taken to pick points which occurred during the distinct 
plateau region of the temperature history since comparisons 
were to be made with steady state calculations. 

The data were subjected to the following criteria, 
and only data which met it were analyzed for mass transfer 
rates. A Nusselt number was computed using film properties 
based upon an average of the measured drop and air tempera- 
tures and the measured test section flow rate. This value 


of the Nusselt number. 


Nu 


calc * 


computed from the Ranz- 


Marshall correlation was compared with an experimentally 


determined Nusselt number, Nu exp , The experimentally 
determined Nusselt number is calculated from Eq. (5.1) 



(5.1) 


where 

X = the heat of vaporization of the liquid 

D = drop diameter 
dm 

— s- = measured mass transfer rate 
d0 

T^ = measured air temperature 

T^ = measured drop temperature 

k = average film thermal conductivity. 

Eq. (5.1) is simply a rearrangement of the energy balance 
for the drop in which the heat transfer coefficient is given 
in terms of the Nusselt number. The criterion is that only 
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data which gives a ratio of Nu ca x G / Nu eX p the range 

0.85-1.15 will be analyzed for mass transfer rates. This 
excluded about one half of the data analyzed. 

Data points which met the criteria are shown on Fig. 

28. They are plotted versus drop Reynolds number, and they 
appear to be randomly scattered about a value of unity. 

Fig. 29 shows the calculated mass transfer rate error 
versus Reynolds number for the data points analyzed. The 
calculated values are the results of steady state film 
theory corrected for high pressure as discussed above in 
this section. The corresponding pressures are indicated 
for each point. All of the points fall within a 30 % error 
band although the average error is - 9 . 1 %. There appears 
to be a systematic trend with Both Reynolds number and pres- 
sure. It is expected that the 50 atm calculated points may 
be high with respect to the 100 atm ones due to the use of 
100 atm interface mole fraction data for both pressures. 

Two explanations of these results are tendered. The 
first is the use of the Ranz-Marshall correlation for drop 
Reynolds numbers above 450. Hugmark (25) recommends a 
slightly different correlation for heat and mass transfer in 
the range of Reynolds number of 450 < Re < 10,000 which 
he claims gives a better fit to the data he analyzed. The 
use of Hugmark* s correlation would lower the three data 
points which have Re > 450 by about 5 % • This would improve 
the consistency of the data somewhat, but it would also 
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Increase the average error from - 9 . 1 % to about -11$. 

The second explanation for the results shown on Fig. 29 
is the influence of drop shape on mass transfer. The low 
Reynolds number drops are slightly more pendant shaped than 
the drops with higher Reynolds numbers. The higher drag 
force associated with higher Reynolds numbers tends to miti- 
gate the deformational effects of gravity since the drag 
force opposes the gravitational force. Using the correla- 
tions of Hsu, Sato and Sage (24) for estimating the influ- 
ence of drop shape on evaporation rates, two trends are 
evident. The first is that the calculated mass flow rates 
at all Reynolds numbers would be 5-10$ higher using the 
Hsu, Sato and Sage results, and the second is that the 
corrections due to drop shape would be somewhat greater at 
the low Reynolds numbers due to greater drop deformation. 

Considering the paucity of data on Fig. 29 and the 
lack of Reynolds number overlap of points at different 
pressures, the explanation given, although plausible, is 
considered provisional. 

Once the mass transfer data was screened by the Nusselt 
number criteria, the effect of pressure was investigated. 

The results of a comparison of calculated mass transfer 
error versus pressure is shown on Fig. 30. This result 
is based upon the 7 data points shown on Fig. 29. There 
are two 100 atm histories, two 50 atm histories, and one 
each at pressures of 10, 5 and 1-1/2 atm abs. The film 
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theory with low pressure assumptions increasingly under- 
estimates mass transfer rates as the pressure is increased. 
At a pressure of 100 atm (reduced pressure of 3*7), low 
pressure film theory underestimates the mass transfer rates 
by 80$. The high pressure film theory greatly improves 
the predicted mass transfer rates at high pressures showing 
an error of only -10# at 100 atm. 

The high pressure film theory shows a small improvement 
with increasing pressure with all the results within 18# 
of the measured values. Over the range of data considered 
the high pressure film theory mass transfer predictions are 
about 10# low on the average. 

Fr-13 Drops Vaporizing in Air 

Approximately 50 individual histories of Fr-13 drops 
vaporizing in air at reduced pressures of 0.75, 1.0, 1.25, 
1.5 and 1.75 were obtained experimentally. Most of the 
data exhibited the classical temperature behavior. The 
drop temperature initially passes from the feed to the 
plateau temperature, resides at the plateau temperature 
until the thermocouple: becomes almost dry, and then rises 
exponentially to the air temperature. At the higher pres- 
sures a few non-classical temperature histories were re- 
corded. All of the data reported here was obtained with 
a chromal-constantan probe thermocouple which had a junction 
diameter of 0.014 in. and 0.003 in. diameter lead wires. 
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Plateau temperatures as a function of air temperature 
are shown on Pig. 31 for a reduced pressure of 0.75 (28.65 
atm abs). Calculated values are given for film theory 
corrected and uncorrected for high pressure. The calcu- 
lated results were obtained by steady state calculations. 

The high pressure film theory contained corrections which 
are summarized in the following list: 

1. Gas phase drop interface mole fractions were 
calculated by the phase equilibrium constant method. 

2. Fr-13 heat of vaporization values were calculated 
by gas and liquid phase partial molal enthalpy 
deviations . 

3. A first order surface regression correction was 
introduced in the form of the air to Fr-13 mass 
flux ratio. 

4. Film thermal conductivity was corrected for high 
pressure by a correlation of Lenoir, Junk and 
Comings (3*0. 

The high pressure film theory underestimates the plateau 
temperature by about 15° F while the uncorrected film 
theory produces only a -10° F error at the higher air 
temperature. 

A series of steady state calculations which show the 
effects of the various corrections were performed for P = 1, 
T^ * 200° F, u =3.7 in. /sec and a diameter of 0.025 in. 

Table V contains the results of these calculations. The 




112 


experimental plateau temperature resulting, from a least 
square fit to the data Is Included for comparison. The 
major contributing correction to increasing the mass trans- 
fer rate is the one for surface regression. The thermal 
conductivity pressure correction is the second most important 
correction from the mass transfer point of view. It increases 
the mass transfer rate by half as much as the surface regres- 
sion correction does. Following the low pressure results, 
the cases listed in the table correspond to the order of the 
list of corrections in the list above. At higher pressures 
the thermal conductivity correction becomes the major fac- 
tor affecting the mass transfer. 

Table V 

Effect of Various High Pressure Corrections 
on Steady State Film Theory Calculated Results 
for Fr-13-Air 

Description 

L. P. film theory 
Film theory (y. ) 

Film theory (y Ao , X) 

Film theory (y Ao > X , N R ) 

H.P. film theory (y Ao ,X,N R ,k' r ) -15 1.000 

Fig. 32 shows similar results for a reduced pressure 
of 1 (28.2 atm abs);as in the case of P^ = 0.75, the low 


Calculated 

Plateau 

Temp . ( °F ) 

-9 

-28 

-22 

-20 


Normalized 
Mass Trans. 
Rate 


0.820 

0.736 

0.601 

0.885 
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Fig. 32. Plateau vs Air Temperature for Fr-13 Drops Vaporizing in Air at 
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pressure film theory gives the best results although typi- 
cal plateau temperature predictions are 5° F low. High 
pressure film theory predicted plateau temperatures are 
about 10° F low. The results for a third calculation with 
all of the high pressure corrections enumerated above except 
for thermodynamic equilibrium at the interface are shown 
on Fig. 32. The phase mole fractions for this case were 
calculated by the phase equilibrium constant method, but 
no air absorption in the liquid phase was allowed. The 
Fr-13 gas phase mole fractions were found to be about 25$ 
higher and the enthalpy of vaporization values 10$ lower 
than the equilibrium values. This assumption gave predicted 
plateau temperatures generally 20° F too low compared with 
the measured values. 

A comparison between high and low pressure theory and 
experimental plateau temperatures is shown on Fig. 33 for 
a reduced pressure of 1.25 (47.8 atm abs). Both calculations 
predict temperatures about 12° F low. However, the low 
pressure theory gives slightly better agreement with the 
measured values at the higher air temperatures. 

Plateau temperature results of measurements and low 
and high pressure film theory calculations at a reduced 
pressure of 1.5 (57.4 atm abs) are compared on Fig. 34. 

There is very little difference between the plateau tempera- 
tures calculated by. either theory. Both theories predict 
temperatures which are about 10° F low. At air temperatures 
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above 2^0° F no plateau temperatures satisfy the high 
pressure film theory. The highest recorded plateau tempera- 
ture at any pressure is plotted on the figure. It corres- 
ponds to a reduced temperature of 0.95 with respect to the 
Fr-13. 

Five data points at a reduced pressure of 1.75 (66.9 
atm abs) were obtained. They are compared with calculated 
plateau temperatures on Fig. 35. There is very little 
difference between the two theoretical predictions except 
that the high pressure theory predicts no solution above 
an air temperature of 160° F. This result is a consequence 
of the property relations; it does not correspond to the 
mixture critical point, and it is believed to be physically 
unrealistic. The calculated values are about 20° F lower 
than the measured plateau temperatures. 

The plateau versus air temperature data were fitted by 
a third order least squares routine and the results for all 
pressures are shown on Fig. 36 . The results appear somewhat 
unsystematic. The curve for a reduced pressure of 1.75 is 
based upon few data, and it may not be representative. How- 
ever, the non-linear behavior of the plateau temperature 
with pressure at low air temperatures is believed to be a 
reflection of the liquid-vapor equilibria characteristics. 

The following figure. Fig. 37, shows this non-linear 
behavior. It is interesting to note that the high pressure 
film theory result has a similar shape to the data which 
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Fig. 37. Summary of Fr-13 Temperature Data, Shown vs 
Pressure for Air Temperature of 150°F 
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were cross-plotted from Fig. 36. The calculated results 
generally predict plateau temperatures 10° F low and become 
more In error at the highest pressure. 

Fig. 38 shows a comparison between the variation of 
plateau temperature with pressure, air temperature fixed 
at 250° F, for calculated and measured results. The high 
pressure theory agrees better with the data at high pres- 
sures, and the low pressure theory gives better agreement 
at low pressures. The calculated results are typically 
10-15° F low over the range of pressures Investigated. 

The Fr-13 data were analyzed for mass transfer rates 
and comparisons were made with calculations. The data 
reduction technique was described in Chapter IV. As In 
the case of the n-heptane drop data, experimentally deter- 
mined Nusselt numbers were compared with ones calculated 
with the Ranz-Marshall correlation at the measured conditions. 
The results for 5 droplet histories at various pressures are 
shown on Fig. 39. The ratio of Nu /Nu are typically 
0.7, and there is no systematic dependence on pressure. 

The range of drop Reynolds number for the points is 275-^25, 
and there is no apparent dependence on it either. 

There are four explanations for the consistently low 


values of Nu , /Nu . The first is that drop motion inclu- 
Cct-i-c exp 

ding internal circulation, oscillation and rotation enhances 


the heat transfer rate with respect :to the rigid sphere 


results of the Ranz-Marshall data. The second is that the 



Plateau Temperature 



0.75 1.0 1.25 1.5 1.75 


Reduced Pressure 


Fig. 38. Summary of Fr-13 Temperature Data, ghown vs 
Pressure for Air Temperature of 250°F 
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Ranz-Marshall correlation corrected for high mass trans- 
fer rates is Inappropriate for vaporizing drops in the near 
critical region irrespective of drop shape and motion ef- 
fects. The third explanation is that the results are pri- 
marily due to the non-spherical shape of the Fr-13 drops. 

As discussed in Chapter IV, in the region of opposed flow 
of combined convection, heat transfer rates are higher 
than those predicted by forced flow for the case of simul- 
taneous heat and mass transfer; this is the fourth explana- 
tion. 

The calculated Nusselt numbers contain a correction 

for high mass transfer rates upon heat transfer as discussed 

in Chapter II. The correction factors varied from 0.65- 

0.8 for the results shown on Fig. 39. The application of 

these calculated correction factors which are based upon 

film averaged properties is a source of error. Application 

of Hsu, Sato and Sage’s results (24) for the influence of 

drop shape on mass transfer indicates that the heat transfer 

rates predicted by the Ranz-Marshall correlation could be 

10-30$ low. However, the Hsu, Sato and Sage correlations 

are for pendant shaped drops, not irregular ones like the 

Fr-13 data, and they are for mass transfer, although the 

heat-mass transfer analogy would support their application 

to heat transfer rates also. Fig. 40 shows the Nu values 

exp 

1/2 1/3 

for the data points normalized by (Re Pr ) plotted 
2 

versus Gr/Re , the combined convection parameter. The 
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results of Narasimhan and Gauvin (38) for aiding and 
opposing flow are shown also on the figure for comparison 
as well as the pure forced convection asymptote of Nu exp / 
(Re 1//2 Pr 1 ^^)=0.68. The Fr-13 drop data fall generally 
between the forced flow limit and that for opposed flow 
supporting the hypothesis that the combined convection 
phenomenon is an important effect. The Fr-13 data may well 
be the result of a starting transient since the droplet 
lifetimes are less than 0.5 sec. Since the results shown 
on Fig. 40 are fairly consistent and there are several 
probable explanations for the low calculated values, the 
data were accepted and mass transfer comparisons were made. 

Figs. 41, 42, 43, 44 and 45 show comparisons between 
quasi-steady film theory calculations and drop histories 
at reduced pressures of 0.75, 1.0, 1.25, 1.5 and 1.75, 
respectively. Drop temperature and mass transferred his- 
tories are compared for each drop. In every case but 
P r = 1, both the low pressure and high pressure film theory 
calculated result is shown. On Fig. 42 the results for 
P r = 1, calculated results for high pressure film theory 
with no absorption of air in the liquid phase, are included. 
In every case the calculated temperatures are lower than 
the measured ones although the calculated values are usually 
within 10° F of the measured. The temperature history calcu- 
lated under the assumptions of high pressure film theory with 
no air absorption in the liquid phase is almost twice as low 
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Fig. H2. History of Fr-13 Drop Vaporizing in Air at P 
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Run 8 (3-14-69) 

Air Temp. = 153°F 
Air Vel. = 2.68 in. /sec 



3.4 3.5 3.6 3.7 3.8 

Time, sec 


Fig. 44. History of Fr-13 Drop Vaporizing in Air at 
!.5 P c 
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Air Temp. = 99 F 
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Fig. 45. History of Fr-13 Drop Vaporizing in Air at 
1.75 P_ 



132 


compared to the measured value as the other calculated 
values are for the same P = 1 case. 

The calculated mass transfer rates are lower than the 
measured rates with the low pressure film theory result 
usually twice as much in error as the high pressure result. 
At reduced pressures from 0.75-1*5 the high pressure film 
theory results are 0 - 30 $ low with the error increasing with 
increasing pressure. The mass transfer rate predicted by 
high pressure theory with no air absorption for the P = 1 
case, shown on Pig. 42, corresponds closely to the measured 
rate although the high pressure theory result assuming equi- 
librium at the interface is only 10$ below the measured 
rate . 

The slopes of the calculated mass transferred curves 
on Pigs. 41-45 were compared with the slopes of the corres- 
ponding measured rates for both high and low pressure film 
theory. These results are shown on Pig. 46. There is a 
distinct trend of underestimation with increasing pressure 
for both theories. Errors of 50$ and 75$ are noted for 
the highest pressure, P^ = 1.75, for high and low pressure 
film theories, respectively. 

The calculation of the critical mixing line or gas 
mixture transition to a critical state was discussed in 
Chapter II. The solid line on Pig. 47 is a portion of the 
locus of critical states calculated for Fr-13-air mixtures. 
As previously discussed, the attainment of a critical state 
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in the gas phase at the drop interface was expected to be 
the initiating condition for non-classical temperature 
history behavior. If this hypothesis is correct, the depar- 
ture from plateau-like behavior should correlate with the 
critical mixing line. This correlation is confirmed by the 
maximum measured drop temperatures exhibiting a distinct 
plateau which are plotted for four pressures on Pig. 47. 

The temperature point corresponding to the three highest 
pressures each resulted from a series of drop temperature 
histories at successively lower air temperatures at a 
fixed pressure. The transition from non-classical to 
classical temperature histories is shown on Fig. 48 for 
reduced pressures of 1.65* 1.75 and 1.78. 

The point on Pig. 47 corresponding to a pressure of 
57.5 atm abs (iP = 1.5) was the highest plateau temperature 
obtained at the pressure. No steadily increasing drop 
temperature histories were recorded since the higher air 
temperatures required would have exceeded the rig maximum 
design limit of 300° P air temperature. 
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VI. CONCLUSIONS 

The suspended droplet vaporization data obtained in 
this investigation are shown in relation to high 
pressure vaporization results of other investigations in 
Pig. 49. The present results are the first vaporization 
data known to the author that were obtained at reduced 
pressures above 0.75 where drop size and temperature his- 
tories were obtained under measured flow conditions. In 
addition, the n-heptane-air vaporization data are the first 
known data at reduced pressures above 2 under any flow 
conditions . 

The conclusions which follow are based upon the 
n-heptane-alr and the Fr-13-air data reported in Chapter V. 
The conditions of the n-heptane data are: air temperature 

100-300° F, air pressure 1.5-100 atm abs, drop Reynolds 
number 150-600, and the initial drop diameter 1275-2025 
microns. The Fr-13 drop data conditions are: air tempera- 

ture 100-300° F, air pressure 29-68 atm abs, drop Reynolds 
number 250-450, and initial drop diameter 980-1200 microns. 
The seven main conclusions are given in the following list. 

1. The ability to suspend vaporizing droplets in 1 g; 
field on a temperature probe at measured liquid 
temperatures up to a reduced value of 0.95 (within 
28° R of the liquid critical temperature) and at an 
air pressure of 56 atm was demonstrated. 



Calculations 



Pig. 49. Summary of High Pressure Drop Vaporization Results 
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2. Under conditions of low ambient gas temperature 
(T /T < 0.8) and high pressure (0 < P/P < 3*7) 

v C 

predictions by film theory corrected and uncorrected 
for high pressure properties both compared favorably 
(within 5° F) with measured plateau temperatures. 

3. Under conditions of low ambient gas temperature 

(T /T < 0.8) and high pressure (0 < P/P < 3.7) 
c c 

predictions by film theory corrected for high 
pressure properties agreed to within + 25 % of 
measured mass transfer rates. Film theory uncor- 
rected for high pressure yielded mass transfer 
results of similar accuracy up to about 5 atm abs 
air pressure (P/P = 0.18); at the highest pres- 

\y 

sure the predicted mass transfer rate was 80% low. 

4. Under conditions of moderate ambient gas tempera- 
ture (1 < T oo /T c < 1.5) and high pressure (0.75 < 

P/P < 1.75) predictions by film theory corrected 

V 

and uncorrected for high pressure properties were 
consistently 10-15° F below the measured plateau 
temperatures. The low pressure theory gave slightly 
better plateau temperature predictions at reduced 
pressures below 1 . 25 , and the results for both 
theories were the same at higher pressures. 

5. Under conditions of moderate ambient gas temperature 

(1 < T /T < 1.5) and high pressure (0.75 < 
c 

P/P < 1.75) predictions by film theory corrected 

c 
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for high pressure properties agreed to within 
35$ of measured mass transfer rates at reduced 
pressures less than 1.5. The predicted mass 
transfer rates were low by about 50$ at a reduced 
pressure of 1.75, and the predicted values exhi- 
bited a systematic trend to increasingly under- 
estimate mass transfer rates as the pressure 
increased. The low pressure film theory mass 
transfer predictions showed similar trends although 
the errors were generally twice as great. 

6. Under conditions of moderate ambient gas tempera- 
ture (1 < T /T < 1.5) and high pressure (0.75 < 
P/P- < 1.75) drop heat transfer Nusselt numbers 

G 

calculated from the Ranz Marshall correlation were 
about 30$ below Nusselt numbers derived from 
the measured vaporization rates and temperature 
differences . 

7. A departure from classical drop temperature his- 
tory behavior was observed at reduced pressures 
greater than 1.5, and this transition from behavior 
exhibiting a plateau to that with steadily increa- 
sing temperature was correlated with the drop 
interface gas mixture critical state. 

In the range of data considered significant . non- 
idealities of thermodynamic and transport properties were 


encountered. There is scant data available and it was 



usually necessary to resort to calculations using general- 
ized equations of state and corresponding state type corre- 
lations of properties. These property non-idealities 
introduced uncertainties into the calculated results and 
indirectly into some of the measured results, e.g., 
experimental Nusselt numbers for heat transfer. 

In Chapter IV the uncertainties in the measured quan- 
tities were estimated as well as their effects upon the 
calculated results. For the Fr-13-alr case considered it 
was estimated that the calculated results should lie within 
13° R for drop plateau temperature and 23 % of mass transfer 
rate of the measured values if the theory itself and pro- 
perty values used within the theory are accurate. All of 
the high pressure film theory calculations agreed with the 
data within these tolerances over the entire range of data 
except with respect to mass transfer rates at reduced pres- 
sures greater than 1.25 and at moderate ambient gas tempera- 
tures. Thus, on this basis the application of film theory 
corrected for high pressure is vindicated for the entire 
range of data of the present investigation shown on Fig. 49 
except for that associated with the Fr-13-air data above 
a reduced pressure of 1 . 25 . 

Corrections to the Nusselt and Sherwood numbers com- 
puted by the Ranz-Marshall correlations for high mass 
transfer rates of 0.65-0.8 were applied to the Fr-13 drop 
calculated histories. These relatively large corrections 
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raise doubts about the application of the Ran'z-Marshall 
results to these data. The corrections to the n-heptane 
calculated results were much closer to unity due to the 
relatively low air temperatures covered. 

There are several recommendations for future work in 
the area of high pressure droplet vaporization. They are 
primarily focused on the approach of suspended droplet 
experiments. The recommendations are as follows: 

1. Additional data with other vaporizing media in 
the high pressure, moderate temperature range is 
desirable. This would allow confirmation of the 
data and film theory calculations developed for 
Fr-13-air. A paraffin fuel vaporizing in nitrogen 
would be a good chpice. Another candidate liquid 
would be CC >2 vaporizing in a nitrogen environment 
since comparisons with Manrique’s detailed calcu- 
lations ( 35 ) would be possible. A paraffin fuel 
or COg vaporizing in high pressure air, although 
more readily obtainable with the present experimental 
setup, would be less desirable but acceptable com- 
binations . 

2. It is desirable to have drop shape measurements 
of unsupported drops moving relative to a gas 
environment at supercritical pressures and tempera- 

’ tures in order to relate the results of suspended 
drop experiments to conditions in a combustion 
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chamber. 

3. The corrections to heat and mass transfer coeffi- 
cients due to high mass transfer rates require 
some attention, perhaps theoretical, since they 
become important at near critical conditions. 

4. High pressure vapor-liquid equilibria data up 
to the mixture critical point is necessary to 
refine any theories and carefully test high pres- 
sure vaporization models with vaporization data. 

5. In the range of low ambient gas temperature 

(T /T < 0.8) and high pressure (0 < P/P < 3.7) 

00 C C 

vaporization data over a greater range of droplet 
Reynolds number at each pressure are desirable to 
check on the Reynolds number dependence of mass 
transfer rate. This would give a greater varia- 
tion in drop shape (independent of pressure) which 
may have a significant effect on the heat and mass 
transfer coefficients. 
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APPENDIX A. PROPERTY DATA 

The property relations for air, water, n-heptane and 
Fr-13 ( chlorotrifluoromethane ) are summarized in this 
appendix. These equations were utilized in the steady- 
state and quasi-steady vaporization calculations. In most 
cases data from the source, indicated by the reference 
number following the property name, were fitted by a cubic 
least squares fit. In some cases property equation was 
taken directly from the data source. The signed number 
following "E" in a coefficient indicates the power of ten 
by which the preceding number should be multiplied. 

Air Properties 

1. Specific heat (22) 

c = (6.2240836 - 9.892961E-03T + 1.2698753E 
P 

-05T 2 - 5.476402E-09T 3 ) 0.068559 
where 504 < T < 756° R and C p (BTU/lb m - 0 R) 

2. Thermal conductivity (22) 

k = ( 6 . 95252 E-O 3 + 2.1973672E-03T - 3.3275553E 
-07T 2 - 6.0211273E-11T 3 ) 0.3225E-06 
where 504 < T. < 756° R and k(BTU/in.-sec-°R) . 
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3. Dynamic viscosity (22) 

y = (-4 . 2220761E-03 + 2 . 57^ 3035E-03T - 1.2586084E 
-06T 2 + 3.621639E-10T 3 ) 0.960E-06 
where 504 < T < 756° R and y ( lb^/in. -sec) 

4. Critical constants (47) 

T c = 132.4° K 
P c = 37-25 atm 

Water Properties 


1. Vapor Pressure (29) 

P ^ - 218.167 x 10" n 
sat 

For T l > 373.16 

n = [3.346313 (647.27 - T^) + 4 . 1413E-02 ( 647 . 27 

- T l ) 2 + 7.515484E-09(647.27 - T L ) 4 + 6.56444E 
-11(647.27 - T L )5] / [ 1.379448lE-02T L 
(647.27 - T l ) + T l J 

For T l £ 373.16 

n =[ 3.2.437814(647.27 - T l ) + 5 . 86826E-0 3(647. 27 

- T l ) 2 + 1.1702379E-08(647.27 - T L ) 4 ] /[2 . 1878462E 
-03 x T l ( 647.27 - T l ) + T l J 

where T^(°K) and P gat (atm) 
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2. Enthalpy of vaporization (29) 

For T l < 690° R 

X = 1.5409974E+03 - 1.5304025 x T L + 1.701721E 
-03 x T l 2 - 1.0238529E-06T L 3 

For T l > 690° R 

X = 408.35902 + 3.1484431T L - 4 . 7l8l765E-03T L 2 
+ 1.9004536E-06T T 3 

li 

where T L (°R) and X(BTU/lb m ) 

3. Vapor specific heat (37) 

C = 0.5799577 - 2.9462E-06T - 1 . 4249794E-06T 2 
+ 1. 81177E-09T 3 

where T(°R) and C p (BTU/lb m -°R) 

4. Vapor thermal conductivity (22) 

k - (-0.15603002 + 2.0814256E-03T + 5.862479E 
-07T 2 - 8.469343E-11T 3 ) 0.212E-06 
where T(°R) and k(BTU/in.-sec-°R) 

5. Vapor dynamic viscosity (22) 

y = (-3.7342072E-02 + 1 . 5360394E-02T + 7.3957384E 
-06T 2 - 3.8423677E-09T 3 ) 0.560E-07 
where T (°R) and y (lb m /in.-sec) 
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6. Binary diffusion coefficient .for water-vapor-air (31) 

„ _ 2.4E-07T 1 * 90 

U AB P 

where T (°R), P (atm) and V ^ (in. 2 /sec) 

7. Liquid density (29) 

P L « (1.5877211E-02 + 1.870l49E-06T L - 1.030489E 

-08T l 2 + l. 393879E-11T l 3 ) 1728 _1 

where T t (°R) and p T (lb /in.^) 

Li j_» m 

N-Heptane Properties 

1. Vapor pressure (45) 

in P sat = 11.94763 - 
where (°R) 

2. Enthalpy of vaporization (45) 

X = 139.9 + 0.181T l - 2. 7875E-04T l 2 

where 500 < T t < 900° R and X(BTU/lb ) 

Li m 

3. Vapor specific heat (45) 

C = 0.119376 + 5.394E-04T 
P 

where 560 < T < 960° R and C p (BTU/lb m -°R) 

4. Vapor thermal conductivity (45) 

k = -1.416E-07 + 6. 2222E-10T 
where 500 < T < 1400° R and k(BTU/in.-sec-°R) 
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5. Vapor dynamic viscosity (45) 

y = 3.92E-08 + 5.88E-10T 
where 510 < T < 1060° R and y (lb m /in.-sec) 

6 . Binary diffusion coefficient for n-heptane-air (45) 

V kB = ("3.131E-03 + 1.2868E-05T + 2.3927E 
-08T 2 )/P 

where 540 < T < 830° R, P (atm), and P^g(in. 2 /sec) 

7. Liquid density (45) 

P L = 3.1662E-02 - 9.5355E-06t l - 6.945E-09T L 2 
where 510 < T r < 650° R and p T (lb /in.^) 

■Li um 

8 . Liquid specific heat (45) 

C T = 0.231 + 5. 62E-04T t 
pL L 

where 500 < T t < 700° R and C T (BTU/lb -°R) 

Li pL m 

9. Critical properties (47) 

T q = 540.3° K 

P„ = 27.0 atm 
c 

Fr-13 Properties 


1. Vapor pressure (16) 

P - ' , 

sat " 14.696 
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for T l < 475° R 

n = 36.7613 - 2623. 988T” 1 - 11.80586 ln(T, )/ 

li 

2.302588 + 5-7l495E-03T L 
for T l > 475° R 

n = 56.34405 - 3351.28lT L _1 - 19.1691 ln(T )/ 
2.302588 + 9.20366E-03 T l 
where (°R) and (atm). 

2. Enthalpy of vaporization ( 16 ) 

for T l < 500° R 

X = -165.41314 + 1.1028303T L - 1 . 336731E-03T L 2 
-1.23949E-07T L 3 

for T l > 500° R 

X = -1251.5327 + 3.766l221T L + 7 . 834936E-04T L 2 
-6.329101E-06T l 3 
where T t (°R) and X(BTU/lb ) 

Ij m 

3. Vapor specific heat (14) 

C = 0.1166976 - 1.0678591E-04T + 4.9459767E 
I? 

-07T 2 - 3.1573555E-10T 3 
where 460 < T. < 760° R and C p (BTU/lb m -°R) 
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4.. Vapor thermal conductivity (14) 

k = [ 0. 005.64 + 1.955E-05CT-460) ] /42300 
where 410 < T < 685° R and K(BTU/in.-sec-°R) 

5. Vapor dynamic viscosity (28) 


U = (-0.017186351 + 1.0192673E-04T - 9.879982E 

-08T 2 + 3 . 7369087E-11T 3 ) 5.60E-05 

where 540 < T < 680° R and u(lb /in. -sec) 

m 

Binary diffusion coefficient for Fr-13-air (4) 

T io (n “ 10) 


where 


n = 7.0013417 + 1.5321532 Log 10 T - 1.5288892 
(Log 10 T) 2 + 0.33510243 (Log 1() T) 3 
and 510 < T < 680° R, P (atm), and V (in. 2 /sec) 


7. Liquid density (16) 

P L - [ 36.07 + 0.01566 (543.59 - T l ) + 3.245E 

-05(5.43.59 - t l ) 2 + 1.11 (543.59 - t l ) 0,5 

+ 6.665 (543.59 - T l ) 1/3 ]/1728 

where 460 < T t . < 540° R and p T (lb /in. 3 ) 

Li ±j in 

8. Liquid specific, heat (15) 


C. = -3.8797818 + 3.4082&99E-02T t - 9.1990079E 
pL n 

-05T t 2 + 8.1380467E-08T 3 

II 1/ 
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where 410 < T t < 535 and (BTU/lb -°R) 
n rL m 

9. Critical constants 

T = 302° K 
c 

P c = 38.2 atm 

N-Heptane-Air Mixture Properties 

1. N-heptane interface gas phase mole fraction (Chapter II) 

y Ao = 9.31717 - 0.00219465 T l - 3 . 7804lE-06T L 2 
+ 8.90531E-09T l 3 

where 540 < < 570° R and P = 1.5 atm abs 

y Ao = exp[-3. 039819 - 0.0565159T L + 1.571044E 
-04T l 3 - 1.034703E-07T L 3 ] 
where 540 < < 720° R and P = 5 atm abs 

y Ao “ exp(4. 86578 - 0.0885904T L + 1. 94609E-04T L 2 
- 1. 16075E-07T l 3 ) 

where 540 < T L < 720° R and P = 10 atm abs 

y Ao ■ -°- 8294.33 + 4.679l68E-03T L - 8. 801203E-06T L 2 
+ 5.59099E-09T L 3 

where 536 < T t < 698 ° R and P = 50, 100 atm abs 

Ju 

2. Enthalpy of vaporization (Chapter II) 

X - 561.331 - 0.603589T L + 1.66334E-03T l 2 
+ 2.49586E-06T 3 
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where 5^0 < T, < 570, P = 1.5 atm abs, and X(BTU/lb ) 

m 

A - (10229.3 - 2.23207T,. - 4 . 12203E-03T t 2 

±J Li 

+ 3.07419E-07) 0.017964 

where 550 < < 650 , P = 5 atm abs, and A(BTU/lb m ) 

X = (8595-66 - 0.654233T l + 2 . 59937E-03T L 2 

- 8.22705E-06T L 3 ) 0.017964 

where 550 < T L < 670 , P = 10 atm abs, and A(BTU/lb m ) 

X = 120.762 - 0.0877l8T L + 2.918424E-04T l 2 

- 3.066951E-07T L 3 

where 536 c < 698 0 R, P = 50 atm abs and A(BTU/lb m ) 

X = 159.652 - 0.5426l4lT L + 1,261746E-03T l 2 
- 9.27936E-07T l 3 

where 536 < < 698 ° R, P = 100 atm abs and A(BTU/lb m ) 

Fr-13-Air Mixture Properties 

1. Fr-13 Interface gas phase mole fraction (Chapter II) 

y Ao = 3.404846 + 1.24854E-02 T l - 1 . 15057E-05T L 2 
+ 6.49599E-09T l 3 

where 420 < T t < ^60° R and P = 28.65 atm abs 

JL 

y Ao = -1.958613 + 0.002193246T L + 1. 04 73HE-05T L 2 
- 8.24939E-09T l 3 

where 430 < T t , < 490° R and P = 38.2 atm abs 
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y Ao = 2 • ^5118 - 9.50644E-03 T l - 1 . 33890E-06T l ‘ 
+ 2.64109E-08T, 3 

JLi 

where 430 < < 490° R and P = 47.8 atm abs 

y Ao = ^ • 53226 - 0.0198138 T l - 1.49072E-05T l 2 
+ 1.823l83E-08T L 3 

where 450 < T L < 515° R and P = 57.3 atm abs 

y Ao = 15.2447 - 0.0413856T L - 4 . 16087E-05T L 2 
+ 1.32523E-07T l 3 

where 450 < T L < 493° R and P = 66.9 atm abs 


2. Enthalpy of vaporization 


X = 16.7298 + 0.0453331T L + 5.20392E-04T l £ 
- 1.10909E-06T L 3 

where 420 < < 460° R and P = 28.65 atm abs 


X = 76.0531 - 0.065l6l7T L - 4 . 23595E-04T l ‘ 
+ 2.033^8E-07T l 3 

where 430 < < 490° R and P = 38.2 atm abs 


X = -104.036 + 0.5l6387T L + 1.91283E-04T l £ 
- 1.43168E-06T l 3 

where. 430 < T,, < 490° R and P = 47.8 atm abs 

. li 


X = -856.017 + 3.305l4T L - 1.30937E-03T l ‘ 
- 3.65246E-06T L 3 

where 450 < < 515° R and P = 57.3 atm abs 
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X = 583.061 - 6.47868T L + 3.69395E-02T L 2 
- 4.33773E-05T l 3 

where 450 < < 493° R and P = 66.9 atm abs. 
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APPENDIX B. DATA SUMMARY 

The water data presented In Chapter V is complete; 
no mass transfer, droplet size or size measurements were 
taken. The Fr-13 histories given in the form of figures 
in Chapter V are complete. However, the n-heptane drop 
mass transfer data in Chapter V are incomplete. A complete 
summary of the n-heptane mass transfer data is contained 
in Table VI. 



Table VI 
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